RENORMALIZATION AND CENTRAL LIMIT 
THEOREM FOR CRITICAL DYNAMICAL SYSTEMS 
WITH WEAK EXTERNAL NOISE 



OLIVER DIAZ-ESPINOSA AND RAFAEL DE LA LEAVE 

Abstract. Wc study of the effect of weak noise on critical one 
dimensional maps; that is, maps with a renormalization theory. 

We establish a one dimensional central limit theorem for weak 
noises and obtain Berry-Esseen estimates for the rate of this con- 
vergence. 

We analyze in detail maps at the accumulation of period dou- 
bling and critical circle maps with golden mean rotation number. 
Using renormalization group methods, we derive scaling relations 
for several features of the effective noise after long times. We use 
these scaling relations to show that the central limit theorem for 
weak noise holds in both examples. 

We note that, for the results presented here, it is essential that 
the maps have parabolic behavior. They are false for hyperbolic 
orbits. 



The goal of this paper is to develop a rigorous renormalization theory 
for weak noise superimposed to one dimensional systems whose orbits 
have some self-similar structure. 

Some examples we consider in detail are period doubling and critical 
circle maps with golden mean rotation number. 

To be more precise, we consider systems of the form 



where / is a map of a one dimensional space (M, or / = [—1, 1]) into 
itself, {^n) is a sequence of real valued independent mean zero random 
variables of comparable sizes, and a > is a small parameter -called 
noise level- that controls the size of the noise. 
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We will study the scaling limit of the effective noise of (11.11) for small 
noise level, and large number of iterations of the system. For / either 
a map at the accumulation of period doubling or a critical circle with 
golden mean rotation and the noises satisfying some mild conditions 
(existence of moments and the like), we will show that the scaling limits 
of the noise resemble a Gaussian in an appropriate sense. That is, if 
we fix a small value of the noise level a and then, look at the effective 
noise after a long time, the distribution of effective noise, normalized 
to have variance one, will have a distribution close to a Gaussian. See 
Theorems 12.11 and 12.21 for precise statements. We will also show that 
statistical properties of the noise, namely Wick-ordered moments or 
cumulants satisfy some scaling properties. Observe that Theorems 12.11 
and 12.21 are similar to the classical central limit and the Berry- Esseen 
theorems, even if the powers which appear as normalization factors are 
different from those in the classical theorems. 

The papers of |CNR81t ISWMSlj considered heuristically a renor- 
malization theory for weak Gaussian noise perturbing one dimensional 
maps at the accumulation of period doubling. The main result in those 
papers was that after appropriately rescaling space and time, the effec- 
tive noise of this renormalized system satisfies some scaling relations. 
The paper |VSK84j developed a rigorous thermodynamic formalism 
for critical maps with period doubling. Among many other results, 
|VSK84] . study the effect of noise on the ergodic theory of these maps 
and showed that for systems at the accumulation of period doubling 
with weak noise, there is a stationary measure depending on the mag- 
nitude of noise that converges to the invariant measure in the attractor. 
A very different rigorous renormalization theory for systems with noise 
is obtained in |CL89j . 

The results we present here can be considered as a rigorous version 
of the theory of |CNR8H ISWM81] . The theory developed here also 
applies to noise of arbitrary shape and shows that the scaling limit is 
Gaussian. The main idea is that one can also renormalize other statis- 
tical properties of the map (cumulants) and, by analyzing the different 
rates of convergence of all these different renormalization operators 
(these are what we call the convexity properties of the spectral radii 
(See Theorem 14. 5p . we obtain that the effective noise is approximately 
Gaussian. The argument we present uses relatively little properties of 
the renormalization operator. Basically we just need that there is some 
convergence to a scaling limit. We also apply similar ideas to the case 
of circle maps. The results there are very similar. 
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Besides the results for the renormahzation, we show that the previous 
results about the renormahzation group give information about the 
behavior of the noise along a whole orbit. 

There is some overlap between the results in this paper and some 
of the results in |VSK84] . The main emphasis on the paper |VSK84] 
is on statistical properties for a fixed level of noise. In this paper, 
we emphasize the behavior on single orbits for longer times but with 
weaker noises. The paper |VSK84j uses mainly the thermodynamic 
formalism and this paper relies mainly on transfer operators. Never- 
theless, there are some relations between the two points of view. Of 
course, the relation between thermodynamics and transfer operators 
goes back to the beginning of thermodynamic formalism. See |May80| . 
In |VSK84t p. 31] the authors introduce one of the cumulant operators 
we use and find a relation between its spectral radius and thermody- 
namic properties. We think it would be possible and interesting to 
develop thermodynamic formalism analogues of the convexity proper- 
ties of transfer operators obtained in Theorem 14.51 It would also be 
interesting to develop thermodynamic formalism analogues of the argu- 
ments developed in Sections 14.61 15.61 which allow to study the behavior 
along a whole orbit from the study of renormahzation operators. 

This paper will be organized as follows. In Section [2] we state a 
general central limit theorem (Theorems 12. ip and a result on Berry- 
Esseen estimates (Theorem l2.2l) for the convergence in the central limit 
theorem for one dimensional dynamical systems with weak noise. The 
main hypothesis of these results is that some combinations of deriva- 
tives grow at certain rate (see (12.91) ). This condition is reminiscent of 
the classical Lindeberg-Lyapunov central limit theorem sums of inde- 
pendent random variables. 

One important class of systems that satisfy the condition (12.91) of 
Theorem 12.11 is that of fixed points of renormahzation operators. Spe- 
cific results for systems at the accumulation of period doubling (The- 
orem 12. 3p and for critical circle maps with golden mean rotation num- 
ber (Theorem [2]3]) will be stated in sections [2^ and [2731 respectively. 

The rest of the paper is devoted to providing proofs of the results 
above. In Section [31 Theorems 12.11 and 12.21 are proved in detail. The 
method of the proof is to use the Lindeberg-Lyapunov central limit 
theorem for a linearized approximation of the effective noise, and then 
to control the error terms. Some examples will also be discussed in 
Section 13. 6[ 

In Section [H we study in detail unimodal maps at the accumulation 
of period-doubhng jFfelTTl [CE801 [CEL80] . and prove Theorem [231 In 
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section we study critical maps of the circle with golden mean rotation 
number and prove Theorem I2.4[ 

The techniques used in Sections [Hand [5] consists on introducing some 
auxiliary linear operators - which we called Lindeberg-Lyapunov op- 
erators (see Sections 14.2. II and I5.3P - that describe the statistical prop- 
erties of the renormalized noise. The crucial part of the argument is to 
show that the spectral radii of these operators satisfy some convexity 
properties (Theorem 14. 5p . so that in the scaling limit, the properties of 
the noise are that of a Gaussian. This implies that the sufficient condi- 
tions of Theorem 12.11 hold for a sequence of times. A separate argument 
developed in Section 14.61 shows that, from the knowledge at these ex- 
ponentially separated times, one can obtain information of the orbit 
starting at zero for all times. From that, one can obtain information 
for orbits starting at any iterate of the critical point. 

2. Statement of results 

Throughout this paper, we will make the following assumptions 
Al. / : M ^ M is a C2 map where M = M, J = [-1, 1] or 
A2. If M = J, we will further assume that there is a number a > 
such that / G C^([-l - a,l + a]) and /([-I - a, 1 + a]) C /. 

For any function / G C(M), we will denote by ||/||co = ^'^PxeM 1/(^)1- 
Let {C,n) be a sequence of independent random variables (defined in 
some probability space (fi, JF, P)), with p> 2 finite moments. We will 
assume that 

A3(p). E[^„] = for all n, and that 

C< ll^nlh < Un\\p < C 

for some constants c, C. Here \\in\\s = (IE[|^nH)^^'*, for any 
s > 0. 

A direct consequence of A3(p) and independence is that 

< n^l^C 

p 

This inequality will be useful to give sufficient conditions on the noise 
level to obtain a central limit theorem (see Theorem 12. II below). 

2.1. General results for one— dimensional dynamical systems 
with external weak noise. Let x G M and a > be fixed. We 
consider the system 

Xm = f{Xm-l) + Xq = X (2.1) 

a will be referred to as the "noise level". 



c < 



max 1^ 
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For M = R or M = T^, define the process Xn{x, a) as the value at 
time n of fl^ . 

For M = I, we define the process a) by (12.11) provided that 
{f{x,_,) + aQ]^,c[-l-a,l + a] 

Observe that 

I J {\f{xj^i) + a^j\ > 1 + a} C la max > a \ := Cn, 

hence, by Chebyshev's inequality we get 

< /[max,<,<„|0|] 
a 

Hence, for a small value of a, the event C„ occurs with low probability. 
Therefore, in order to define x„(a;, cr), it suffices to condition on the 
event 



^\Cn = {cr max < a 

{ l<i<n 

In the scaling limit, we will consider a sequence of noise levels {cr„} 
converging to so that f2 \ C„ is close to the whole space fl. Notice 
that if (^j)jeN is supported on a compact interval, then by taking cr^ 
small enough, the events C„ will be empty sets. 

2.1.1. General central limit theorem for one dimensional dynamical sys- 
tems with random weak noise. Our first result is a general central limit 
theorem for one-dimensional dynamical systems with weak random 
noise (Theorem l2.1l) . This result is based upon the classical Lindeberg- 
Lyapunov central limit theorem for sums of random variables }Bil68l p. 
44]. This is reflected in the sufficient condition (12. 9p in the statement 
of Theorem 12.11 

We introduce the following notation 

Definition 2.1. Let / be a map on M satisfying Al, A2, and let (.^„) 
be a sequence of random variables with p > 2 finite moments, that 
satisfies A3(p). 

1. The Lyapunov functions As{x,n), for s > 0, and A are defined 
by 

n 

As{x,n) = 5^|(r-0'of (x)[ (2.2) 
i=i 

i 

A{x,n) = maxE|(/"'T°/'(^)| (2.3) 
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When needed, we will use the notation A-^ and A{ to emphasize 
the dependence on /. 
2. Let X E M and o" > be fixed. The linearized effective noise is 
defined as 

n 

Ln{x) = Y,{r'')'or{x)^, (2.4) 

Remark 2.1. It is very important to observe that for each s > 0, the 
Lyapunov function As{x,n) satisfies 

A,(a;, n + m) = \ (D' o /"(x) 1' A,(x, n) + A,(r(x), m) (2.5) 

We will use (12. 5p in the study of central limit theorems for systems 
near the accumulation of period doubling, Section HI and for critical 
maps of the circle. Section [51 

Remark 2.2. The Lyapunov functions (12. 2p are used to estimate the 
sums of the moments or order s of the terms in (12.41) . In particular, 
notice that 

var[L„(x)] = ^ ((r^^)' o f (x)) E[e|] (2.6) 
i=i 

The assumption A3(p) and Holder's inequality imply that there are 
constants c, C such that for any < s < p, 

n 

cKsM < ^ |(r-^yof (a;)[E[|0l1 < C^M (2.7) 

In particular, if A3(p = 2) implies 

cA2(x,n) < var[L„(x)] < CA2(x,n) (2.8) 

The main result that we obtain for orbits of a dynamical system 
is that if a deterministic condition on the orbit (expressed in terms 
of Lyapunov functions) holds, then, the noise perturbing this orbit 
satisfies a central limit theorem. 

Theorem 2.1. Let f , M be a function satisfying Al and A2, and let 
(^„) he a sequence of independent random variables, with p > 2 finite 
moments, that satisfies A3(p). Suppose that for some a; G M there is 
an increasing sequence of positive integers Uk such that 

lim ' = 2.9 

'^^^iA2ix,n,)f' 

Let {ak)k be a sequence of positive numbers. Assume furthermore 
either of the two conditions 



CLT FOR DYNAMICAL SYSTEMS WITH WEAK RANDOM NOISE 



7 



HI The noise satisfies A3{p) with p > 2. 
And the sequence (cjfc) satisfies: 

j.^ |inico||maxi<j<„, \^j\\\l(A{x,nk)fak ^ ^ ^g) 

H2 The noise satisfies A3{p) with p > 4:. 
And the sequence (ak) satisfies: 

hm = (2.11) 

Then, there exists a sequence of events E J-" such that 

Ml. limfc^ooP[5fc] = 1 

M2. The processes defined by 

= (-M-f'^if (2.12) 

Wn,{x,ak) = ^ ' ^ ^ (2.13) 

Vvar[(x„,(x,afc) - /'^fe(x))lBj 

converge in distribution to a standard Gaussian as k ^ oo. 

Furthermore, if the sequence ^„ is supported on a compact set then, 
we can choose = fl for all k. 

The sets Q\Bk, which we call outliers, are events where large fluc- 
tuations of noise occur. We will refer to 



ak) - f-Ax))lB, (2.14) 



as the effective noise. 

Condition (12. 9p in Theorem 12.11 is closely related to the Lyapunov 
condition of the classical limit theorem for sums of independent random 
variables applied to the linearized effective noise Ln^x) deflned by 

We will show in Section [373] that if the sequence of noise levels (cTfe) 
satisfy (12.101) then, the events where large fluctuations of noise occur 
(outliers) have low probability. As a consequence, we will have that 
the linearized effective noise (12.41) is, with large probability, a very 
good approximation to the effective noise deflned by (I2.14p . 

Remark 2.3. Note that there are two variants of the results (12.121) 
and (I2.13P in Theorem 12.11 as well as two variants on the hypothesis. 
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The difference in the conclusions is that the effective noise is nor- 
mahzed in two - in principle different - ways. The version (12.121) nor- 
malizes the effective noise by its variance and the version normalizes 
by the variance of the linear approximation (see Remark 12.21 

It could, in principle happen that the difference between the linear 
approximation and true process converged to zero in probability but 
that had a significance contribution to the variance. This pathology can 
be excluded by assuming that the noise has sufficient moments and that 
the noise is weak enough. These two hypothesis can be traded off. In 
the first version of the hypothesis [HI], we use only p > 2 moments and 
a somewhat stronger smallness conditions in the size of noise (12.101) . In 
the second version of the hypothesis [H2], we assume p > 4 moments, 
but the smallness conditions in the noise are weaker. 

The subsequent Theorems 12.21 will be true under either of the hy- 
pothesis. 

Remark 2.4. As we will see later, in the proof, the conditions on 0"^ 
are just upper bounds. If we consider two sequences < ak and 0"^ 
satisfies either of (12.101) or (12.111) . then dk satisfies the same conditions. 

Furthermore, the bounds that we obtain for the proximity Wn,,{x, dk) 
to the standard Gaussian are smaller than the bounds that we obtain 
on the proximity of w„j,(a;, ak) to the standard Gaussian. 

2.1.2. Berry-Esseen estimates. The relevance of the outliers will be- 
come more evident in our second result. Theorem 12.21 below, which 
provides the rate of convergence to Gaussian in Theorem 12. 1[ 

We will use $(z) to denote the distribution function of the standard 
Gaussian measure on (M, i3(]R)), that is 

$(z) = ^ r e-^'/^dt 

J — oo 

Theorem 2.2. Let f , {C,n)n be as in Theorem \2.1\ and let s = mm{p, 3). 
Assume that condition (12.91) holds at some x E M If ak is a sequence 
of positive numbers such that 



{A{x,nk)f 
v/var[L„,(x)] 

then, we have that 



Co 



max I C, 



i<j<nk 




As{x, Uk) 
iA2ix,nk)f' 



2 

(2.15) 



sup |PK,(a:)lB, <z]- ^z)\ < A^At^l^ (2.16) 



\A2ix, nk)y/^ 



where A > depends only on x. 
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Remark 2.5. When M = M, or T-*^ and f{x) = x + c for some con- 
stant c, Theorems 12.11 and 12.21 coincide with the classical central limit 
and Berry-Esseen theorems for sums of independent random variables. 
Since /" = 0, the result holds regardless of what values ak takes, and 
the outliers are empty sets. 

Remark 2.6. In section [3] we will consider the map f{x) = 2x. This 
system does not satisfy (12. 9p . and indeed the conclusion of Theorem 12. II 
fails. Systems with enough hyperbolicity satisfy other types of central 
limit theorems for weak noise |GK97] . or even in the absence noise 
[Liv96j . [FMNT05] . Those results are very different from the ones we 
consider in this paper. 

2.1.3. Sketch of the proof of Theorems \2.1\ and \2.^ The proof of these 
results is obtained in Section [3] by showing that: 

1) The linear approximation to the process satisfies a central limit 
theorem (or a Berry-Esseen theorem) 

2) Under smallness conditions on the noise level (see (I2.10p ). the 
linear approximation is much larger than the Taylor reminder, 
so that we can transfer the Gaussian behavior from one to the 
other. 

The main source of difficulties in the proof are situations when the 
noise is much larger than expected from the statistical properties of 
the linear approximation (outliers). These events, of course have small 
probability and, therefore, do not affect the convergence in probability. 
However, it could happen in principle that they change the variance. 
We will see in Section 13.3.31 that the variance of the effective noise 
and that of its linear approximation are asymptotically equal. As a 
consequence, we will have the effective noise normalized by its variance 
converges in distribution to the standard Gaussian. 

The existence of moments of high order will provide an improvement 
on the choice of the noise level. In particular, we will see that if the 
noise has compact support, large fluctuations of the effective noise never 
occur; that is, the outliers are empty sets. 

The procedure of cutting off outliers in Theorem 12.11 is similar to 
the process of elimination of "Large fields" that occurs in the rigorous 
study of renormalization group in |GK85l IGKK87j . 

In this paper, we consider two examples of maps that have a renor- 
malization theory, and for which Theorems 12.11 and 12.21 apply. Namely, 
systems at the accumulation of period doubling and critical circle maps 
with golden mean rotation number. 
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2.2. Results for systems at the accumulation of period dou- 
bling. In section 14.11 we consider systems of the form (11. ip where / is 
a 2A;-order analytic unimodal map of the interval I onto itself. That is 
/(O) = 1, /(J)(0) = 0, for 1 < 2fc - 1, /(2*^)(0) ^ 0, and xf'{x) < for 
X 7^ 0. Here, /'•■^•' denotes the j-th order derivative of /. 

The period doubling renormalization group operator T acting on the 
space of unimodal maps is defined by 

T/(x) = f{Xfx)/Xf 

where \f = /(I). We refer to Section (14.11) for a precise definition of 
unimodal maps and the period-doubling renormalization operator. 

For each k, there is a set of analytic functions Wsigk), such for maps 
/ e yVs{gk), we have that T"/ converges to a universal function 
which is a fixed point of T (see jFfelTTl ILan82l [TC781 |Eps86i ISuM 
IMar98l[dM^^|jS02] l 

We will show in Section H] that 

Theorem 2.3. Let f G Wsigk) then, 

Fl. For any x = f\0), I G N and any p > 2, 

hm ^^^^'^^ = 

n^oo {A2(x,n)}P/2 

Let {C,n) be a sequence of independent random variables which have 
p> 2 finite moments, and that satisfies A3(p). 

F2. Let Wn{x, a) be 
Wn{x,a} - 



Then, there is a constant 7 > 0, such that if 
lim cr^n''''''^ = 

n 

then, for each x = f'' (0) , I E N there are events {i?„(x)}„ C 
with 

limP[5„(x)] = 1 

n 

such that Wn{x,an)'i-Bn(x) and hence Wn{x,an) converge in dis- 
tribution to the standard Gaussian. 
F3. Furthermore, there are constants a > and 6 > depending 
on p such that if cfn < , then 

sup \F[wn{x,an)lB^{x) < z] - <^{z)\ < C^u-'^ (2.17) 

zeK 
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Some explicit values for 7, a, and 6 that follow from the arguments 
of the proof are given in Section 14.61 

In |DEdlL06] . empirical values for 7 are obtained for the quadratic 
Feigenbaum fixed point. The values obtained there suggest that the 
sequence of level of noise decays as a power 7^, of the number of 
iterations, that is not very different from the one we obtain in Section 
SH 

2.3. Results for critical circle maps. In Section [5] we study critical 
maps of the circle with golden mean rotation number [Lan84l Eps89 



IdFdMQQ] . That is, we consider strictly increasing analytic maps / such 
that fix + 1) = fix) + 1, /(^^(O) = for J = 1 . . . , 2A;, f(^^+^\x) ^ 0, 
and lim„(/"(x) - x)/n = (^5 - l)/2 = (3. 

2.3.1. Central limit theorem for Fibonacci times. Recall that the se- 
quence of Fibonacci numbers iQn), is defined by Qo = 1 = Qi, Qn+i = 
Qn + Qn~i- Any integer n admits a unique Fibonacci decomposition 

where mo > . . . > mr„ > are non-consecutive integers, (i. e. rrii > 
rrii-i + 2). Notice that r„ + 1 is the number of terms in the Fibonacci 
representation of n, and that r„ < mo < [log^-i n]+ (Here [ ] stands 
for the integer part function). 

Theorem 2.4. Let f be a critical circle map. If {uk} is a increasing 
sequence of integers such that 

lim — = (2.18) 

k-^oo log^-i rifc 

then, 

CI. For all x = f\0), I G N, and any p> 2, 

hm ^^^^'^^^ = 
fc-^ {A2(a;,nfc)}p/2 

Let i^n) be a sequence of random independent variables that have p > 2 
moments, and that satisfies A3(p). 

C2. Let Wnix,a) be the process defined by 

(x„(x,(t) - /"(x)) 

Wnix,n) = = — 

(jA/var[L„(x)] 

Then, there is a constant 7 > such that if 

limakuV'^ = 
k 
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there are events {Bk{x)} C with 
\imF[Bk{x)] = 1 

k 

such that Wn,^{x,ak)'i-Bk(x) andhenceWn^^{x,ak) converge in dis- 
tribution to the standard Gaussian. 
C3. Furthermore, there are constants r > and v > depending 
on p such that if cxk < "^fc^^; then 

sup \F[wn^{x,ak)lB,{x) <z]- (^{z)\ < D,^n-" 

2.3.2. Central limit theorem along the whole sequence of times. In Sec- 
tion 15.6.31 it is shown that the Lyapunov condition fl2.9p . for orbits 
starting in {/'(O) : / G N}, holds along the whole sequence of times 
provided that some numerical condition fl5.44p is satisfied (see Propo- 
sition [5]6] and Theorem 15. 7p . The technical condition fl5.44p is some 
relation between properties of fixed points of a renormalization oper- 
ators and the spectral radii of some auxiliary operators. It could be 
verified by some finite computation or implied by some monotonicity 
properties of the fixed point theorem, 

Under the technical condition (15 ■441) . the central limit theorem holds 
along the whole sequence of times. The proof presented gives results 
that do not depend on the condition (15.441) . Namely, we show that 
there as a central limit theorem along sequences of numbers which can 
be expressed as sum of sufficiently "few" Fibonacci numbers in terms 
of their size. The hypothesis (I5.44p implies that all the numbers satisfy 
this property. 

Theorem 2.5. Let f be a critical circle map. Under the numerical 
condition (15.440 (see Section [5. 6. 3\) 

C4. Ifx = /'(O), I eN andp> 2, then 

hm ^^^^'^^ = 

n^oo {A2{X, n)}P/^ 

Let (^„) be a sequence of random variables which have p > 2 moments, 
and that satisfies A3. Then, 

C5. There is a constant 6 > such that if 

lim cr^n''''''^ = 

k 

there are events {Bn{x)} C J-" with 
\imF[Bn{x)] = 1 

n 

such that Wn{x,ak)lB„(x) and hence Wn{x,an) converge in dis- 
tribution to the standard Gaussian. 
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C3. Furthermore, there are constants r > and v > depending 
on p such that if cr„ < , then 

sup|P['«;„(a;,(T„)lB„(^) < z] ~ $(z)| < D^n-"" 



2.3.3. Discussion of the results for critical maps (Theorems \2.3[ 2.4 



and \2.5\} . The proof of Theorems 12.31 and 12.41 are obtained using renor- 
mahzation methods. 

Roughly, renormahzation gives us control of the effects of noise on 
small scales around the critical point for a fixed increasing sequence 
of times (powers of 2 in the case of period doubling and Fibonacci 
numbers in the case of circle maps with golden mean rotation). This 
gives us, rather straightforwardly, a central limit theorem when the 
orbit of zero is observed along these sequences of times. 

To obtain a central limit theorem along the sequence of all times, 
we use the fact that an arbitrary number can be written as sum of 
these good numbers. We argue by approximation. We observe that 
the sequence of times accessible to renormahzation is also the sequence 
of times at which the orbit of zero comes close to the origin. Hence, 
we can write the orbit of zero as a sum of approximate Gaussians. 

The argument we present has some delicate steps. We need to bal- 
ance how close is the approximation to the Gaussian (how fast is the 
convergence to the CLT) with how fast is the recurrence at the indi- 
cated times. 

In our approach, to prove a CLT along all times, we have to compute 
and compare the two effects. This comparison depends on quantitative 
properties of the fixed point of the renormahzation group and some of 
the auxiliary operators. 

In the period doubling case, the properties required by our approach 
can be established and proved by conceptual methods (convexity and 
the like) from the properties of the fixed point. 

In the case of circle maps however, our methods require a property 
(see f l5.44l) ) which seems to be true numerically, but which we do not 
know how to verify using only analytical methods. 

The analysis presented above raises the possibility that, for some 
systems, the weak noise limit could have a CLT along some sequences 
but not along other ones. Of course, it is possible that there are other 
methods of proof that do not require such comparisons. We think that 
it would be interesting either to develop a proof that does not require 
these conditions or to present an example of a system whose weak noise 
limit converges to Gaussian along a sequence of times but not others. 
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3. Proof of Theorem 12.11 and Theorem 12.21 

In this section, we prove the two general theorems about one dimen- 
sional dynamical systems with weak random noise, namely Theorem 
12.11 (a central limit theorem) and Theorem 12.21 (a Berry-Esseen theo- 
rem). 

First, we consider a linear approximation of the system (11.11) and 
show that it satisfies a central limit theorem and a Berry-Esseen the- 
orem, see section 13. 1[ Then, in section 13.21 we make a comparison 
between the linear approximation process and cr), see Lemma 

13. 2[ The proofs of Theorems 12.11 and 12.21 are given in Sections 13.31 and 
13.41 respectively. We show that for cr small enough, see (12.101) . the lin- 
ear approximation process is a good approximation to a), except 
perhaps in sets of decreasing probability, which we call outliers. Since 
these sets have probability going to zero, they do not affect the con- 
vergence in probability. Showing that these outliers do not affect the 
variance requires some extra arguments, which we present later. 

3.1. Linear approximation of the effective noise. For x G M and 

cr > fixed, using Taylor expansion, we decompose the process Xn{x, a) 
as 

Xn{x, a) = /"(x) + crL„(x) + a^Qn{x, cr), (3.1) 

where the linear term L„ is the sum of independent random variables 
defined by (12. 4p . The linear approximation process yk{x,a) 

yn{x,a) = r{x)+aLn{x) (3.2) 

satisfies the following central limit theorem. 

Lemma 3.1. Let f he a function satisfying Al and A2, and {S,n) be a 
sequence of independent random variables with p > 2 moments, satis- 
fying A3. If condition (12.91) holds for some point x E M then, 

lnSx)^^h^p= (3.3) 
A/var[L„(x)J 

converges in distribution to the standard Gaussian as A; —> oo. More- 
over, there is a universal constant C such that 

\n.M < .] - < c ^/;'°''''''Xt, , (3.4) 

zdR {A2{x,nk)) 
Proof. From (12.71) and (12. 8p . there is c > such that 



y/A2{x, Uk) v^var[L„^(a;)] y/A2{x, Uk) 
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Hence, the convergence of to the standard Gaussian follows from 
the classical Lindeberg-Lyapunov central limit theorem [BilGSj p. 44], 
since condition (12.91) is equivalent in this case to the Lyapunov condition 
for sums of independent random variables. 

The second assertion, (13.41) . follows from Berry-Esseen's Theorem 
for sums of independent random variables |Pet75l p. 115] □ 

Remark 3.1. Notice that the linear dependence of yn{x,cr) on the 
random variables ^„ implies that the convergence in Lemma 13.11 is in- 
dependent of the noise level a. The size of the noise a will be important 
in the control of the non-linear term a'^Qn{x, a). 

3.2. Nonlinear theory. In this section, we prove a result (Lemma 
13.21) that will help us make a comparison between x„(x, a) and yn{x, a). 
This result is analog to a well known result on variational equations 
for ODE's [Har82j . The method we use in the proof of Lemma [3.21 is 
similar to the proof of the Shadowing Lemma in [Shu87l IKat72l [SziiS 1] . 

In the rest of this section, we will use the norm ||x|| := maxi<j<m \xj\ 
for vectors x in Euclidean space M*", and the corresponding induced 
norms for linear and bilinear operators. 

Lemma 3.2. Suppose f G C^iM), xq G M, and let A e M^+^ with 
Ao = 0. Consider the sequences x, x and c in M^+^ defined by 

a) Xq = Xq, Xi+i = f{Xi), 

b) Xi+i = f{xi) + Ai+i, 

c) Q+i = f {xi)ci + Ai+i with Co = 0, 

fori = 0,...,N -1. 

Assume that 

||A||||r||co{A(xo,iV)r <^ (3.5) 

Then, 

||x + c-x||<||cfA(xo,iV)||r||co 

<||Af(A(xo,iV))3||r||co 

Proof. Fixing a parameterization on M, we can assume without loss 
of generality that M = M. We will prove (13.61) by showing that x is a 
fixed point of the contractive function ip defined in (I3.16p . 
Notice that 

i 

c^ = J2{M{x,)A, 

j=0 
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If A = {ttij) is the lower triangular matrix defined by 

(/'"■')' (xj) for < j < i < 
otherwise 



we have that 

c = AA (3.7) 

||A|| = A(xo,iV) 

Let us define a function r : M^+^ x M^+i — > R^+i by 

r(z, a) = (xo, f{zo) + ai, . . . , f{zN-i) + 

Observe that 

a) r(x, 0) = X 

h) t{z, A) = z if and only if z = x 

c) For all (z, a) G x R^+i 

Dr(z,«) = Dr(z,0) 

Dnr(z,a) = Ait(z,0) 

The following identities will be useful 

L>ir(x,0)c + D2i-(x,0)A = c (3.8) 

(A(r(x,0))-/)-i = A (3.9) 

A direct computation shows that for any point (z, a) and vectors [h, k], 
[h, ~k] G M^+i X R^+i 

Z}V(z,a)([/i,A;],[/i,^]) = (0, /iq/" (^o)/io, • • • , /^TV-i/" (^^-i)/i^-i)^ 
Therefore 

sup ||DMz,«)|| = sup ||Dnr(z, 0)11 < ||r||co (3.10) 

(z,a) z 

We define an auxiliary function Af : M^+i x M^+i — > M^+i by 
J\f{z, a) = -A(r(z, «) - z) + z 

Using (a)-(c) we have 

d) A/'(z, 0) = z if and only if z = x 

e) A/'(z, A) = z if and only if z = x 

/) DiAf{z,a) = DiAf{z,0) for all (z,a). 

It follows from ([SI]) 

M{z,a) -N'{z,0) = -A[r(z,a) -r(z,0)] (3.11) 
Di^{z,a) = -A(Dir(z,0) -/) + / 

(3.12) 

= A[Z}ir(x,0)-Dir(z,0)] 
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D2Af{z,a)k = -A[0,ki,...,kN] (3.13) 
From (13.81) we have 

r(x + c, A) = X + c + ^DuTii, A)(c, c), (3.14) 

where (^, A) is a point on the hnear segment between (x + c, A) and 
(x,0). Thus, by ([321), ISTTOD . and (KWf we have 

||Ar(x + c,A)-(x + c)|| < i||A||||c|nir||co 

< i||AfiiAiniriico 

Let -B(x; r) be the closed ball in centered at x with radius r = 

2||A||||A||. Then, by flXT^ . we have that for all {h,a) G M^+^ x R^+i 
and all z G -B(x;r), there is t] in the linear segment joining z and x 
such that 

Dl^^{z, a)h = -ADur{r], 0){h, z - x) 
Therefore, provided that condition (13. 5p holds, we have that 

||Z}iAr(z,a)||<2||A|n|A||||r||co<^ (3-15) 

for all z G 5(x; r). 

On the other hand (13.111) and (I3.15P imply that 

||A^(z,A)-x|| < ||Ar(z,A)-AA(x,A)|| + ||Ar(x,A)-Ar(x,0)|| 

1,, -n n . MM . n 

< -||2; — x|| + II A|| II A|| < r 

It follows that the function 

tlj{z) = Af{z, A) (3.16) 

is a contraction by a factor 1/2 of the closed ball B{^\r) into itself. 
Moreover, x is the unique fixed point. 
Inequality (13.61) follows from 

||x + c-x|| < ||x + c-Ar(x + c, A)|| + ||A/'(x, A) -A/'(x + c, A)|| 

< i||A||||cf||f'||c;„ + i||x + C-x|| 

□ 
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3.3. Proof of central limit theorem (Theorem 12. The conclu- 
sion of Theorem 12. II will follow by combining Lemmas 13.11 and 13 . 2[ The 

key argument is to define events (outliers) where the linear approxima- 
tion process is very different from x„(x, o"), and show that they have 
small probability. 

Notice that if ||/"||co = 0, then Theorem 12.11 coincides with the 
Lindeberg-Lyapunov central limit theorem. Therefore, we will assume 
that liriico >0. 

3.3.1. Outliers. For each k and j = 1, . . . , n/c, let Aj = a^^j, and let A, 
X, X and c be as in Lemma [221 Then, by (12.51) . (13. 2p and the definition 
of the linear approximation process, we have that for j = 1, . . . , 



yj{x,(Tk) 

(JkLj{x) 

Xj{x,(Tk) 



Xj + Cj 



For any sequence {a^) of noise levels decreasing to 0, Lemma \^72\ implies 
that in the event 



||/"||co(Tfc(A(a;,nfc))2 max \^j\ < ^ 



(3.17) 



the linear approximation process yn^{x,ak) (13.21) is close to x„^(a;, 0"^). 
Thus, we will restrict the process x„^(x, 0"^) to events & J-" such that 



Bu 



max \C,j\ < 



C B, 



(3.18) 



for some appropriate sequence a^. We will choose (a^) so that 

lim F[n \Bk]^0 (3.19) 



For this purpose, it will be enough to define (a^) by 

1 



where j3 and a are chosen so that (I3.18P and (13.191) hold. This means 
that 



ak{A{x,nk)y <a^{A{x,nk)y 

l/a 



lim 

fe^oo 



max l^j 



l<i<nfc 



(A(x,nfc))^/Vfc = 



(3.21) 
(3.22) 
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Remark 3.2. Recall that if M = /, we define the process Xnf,{x,ak) 
conditioned on the event 



\ Cn^ = < ak max |6| < a 

y i<j<nk 

Chebyshev's inequality implies that 

Notice that (I322D, and the fact that A(a;, Uk) > 1 imply that P[C„J 
as — + oo. 

We refer to {Q \ Bk}k ({^ \ (-Bfe fl Ck}k if M = J) as the sequence 
of outliers. 

3.3.2. Estimates on the effective noise. In this section, we will show 
that the outliers have small probability. This means that linearized 
effective noise (12.41) is a good approximation of the effective noise (12.141) 
with high probability. As a consequence, we will have that the effective 
noise scaled by the standard deviation of the hnearized effective noise 
approaches a Gaussian. 

Recall from (13.11) that the effective noise is decomposed as 

To control the effect of nonlinear terms in the effective noise, it will be 
enough to require that the variance of the noise of crlQn^^x, ak) is small 
compared to the variance of the linearized effective noise akLn^{x). In 
terms of the scaled processes w„j,(x,crfe) and luki^) defined by (I2.12p 
and (13.31) . this requirement is equivalent to 

lim \\{wn^{x,ak) - lndx))lBk\\2 = (3.23) 

fc— >oo 

By (I3.18P and Lemma [3.21 we have that 
|Qnfc(a;,cTfc)|lij 

Then, by ( Km we get 

(\(r 77, 1 "1 3-2/3 l-2a 

||K,(x,a,) - ln,{x))lsj2 < C^^^7==4^ (3.24) 

^y A2{X, Uk) 

for some C > 0. To obtain (I3.23p . it suffices to require that 
, (A(a;,nfc))3-2/3a^-2" 

lim ^ ^ ' /^^ — ^ = (3.25) 

VA2(x,nfc) 

We have the following Theorem 



\Qn,{x,ak)\lB, < ||/"||c;o(A(x,nfc))3 max l^^fls, 

i<j<"fc 
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Theorem 3.3. Let f and be as in Theorem \2.1\ Assume that 
the Lyapunov condition fl2.9p holds and that (ak) satisfies fl2.10p . // 
a'-^ = 2 = (3, then Wn^{x,ak) converges in distribution to the standard 
Gaussian. 

Proof. Notice that the Lyapunov condition (12. 9p imphes that 

hm A2(x, oo 

A:— >oo 

Hence, by fl3.2ip and (13.250 . it wiU be enough to consider < a < 1/2 
and (3>2. 

In particular, if = 2 = [3 and (cTfc) satisfies (l2.1Up . then 1^^ 
and (w„j.(x, (jfc) — converge in Lp/2(IP) and hence in prob- 

ability to 1 and respectively. By Lemma 13.31 and the "converging 
together" Lemma |Dur05l p. 89], it follows that w„j,(a;, 0"^) converges 
in distribution to the standard Gaussian. □ 

The result of Theorem 13.31 is improved by considering the weaker 
condition (12. lip for {ak) and smaller outliers Q\Bk. 
Notice that by Chebyshev's inequality 

m\Bk] < (4||/||co(A(x,n,))2||maxi<,<„je„|||pCrfc)'' (3.26) 

On the other hand, we have that 

\/var[L„,(x)] 

The following result is a direct consequence of the ^^convergence to- 
gether" Lemma, I^M>, and IK27} . 

Corollary 3.4. Let f and (^„) be as in Theorem \2.1\ Assume that the 
Lyapunov condition (12. 9p holds and that (ak) is a sequence of positive 
numbers that satisfy (12. lip . Then, the process Wn^i^x, crfc)!^^ and hence 
Wn^{x,ak) converge in distribution to the standard Gaussian. 

3.3.3. Gomparison between the variance of the effective noise and the 
variance of the linear approximation. Observed that the scaling given 
in (I3.24P uses the variance of a random variable L„^(x) which is defined 
in the whole space fl. Since the outliers occur with very low probability, 
it is more natural to restrict all the quantities to the complement of 
the outliers, that is to the events Bk. 

If (cTfe) satisfies (I2.10p . we will show that the variance of the effective 
noise and that if its linear approximation are asymptotically equal. By 
the converging together Lemma, we will have as a consequence a central 
limit theorem for the effective noise scaled by its variance. 



^ ^ II maxi<j<„^ WljAjx, nk))^ak ^7) 
„/2 ~ ^/A2ix, Uk) 
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Lemma 3.5. Let f and (^„) be as in Theorem \2.1\ Assume that the 
Lyapunov condition (I2.9p holds. For any sequence (dk) of positive num- 
bers such that 

lim4-(-^)i5(^^=0, (3.28) 
define the event Dk = [niaxi<j<„j. < dk]. Then, 

lim = 1 (3.29) 

Proof. Notice that 
where 

First, we estimate cov[,^jl£)j., ^jId^] when j ^ i. 

Since IE[^j] = for all j, from the independence of (,^„.) we get 

m,iD,]\ = m.m < dk}]\\[nu\ < dk] 
= m,m>dk}]\iinu<dk] 

By Holder's and Chebyshev's inequality we get 

mM < uxdl'' Umn.\ < 4] (3.30) 

Notice that fl3.28p implies that limfcc/fcoo and that limj!cP[Dfc] 1- 
Consequently, the product of probabilities on the right hand side of 
fl3.30p is close to (and smaller than) 1. Furthermore, A{p) implies that 
(^„) is bounded in Lp(P). Hence, 

mM<Cdl-^ (3.31) 
Similarly, the independence of (^„) implies 

mj^dD,]\ = m,m > dk}]mm > 4}]i n ^ti^^i ^ 
= m,m > dk}mim > 4}]i n ^ti^-i ^ ^'^i 
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Using Holder's and Chebyshev's inequalities we get 

\m^^^D,]\ < m'Mm'^'-'^ n ^ d>^] (3-32) 

rn^j,i 

The product of probabilities on the right hand side of (13.321) is close to 
(and smaller than) 1. Therefore, using A3(p), we have: 

m,^^^D,]\<Cdf^'~''' (3.33) 

Combining (13.311) and (13.33^ we have, for j ^ i, that 

\cov[^AD,,^dD,]\<Cd,'^'-'^ (3.34) 

When j = i, we obtain 

|var[e,lz.J-var[e]| = m^ls^] + mAo,]]' 

< C{¥[Dl]^P'^^^P + df^''-^^) (3.35) 

Recall that A2{x,nk) and var[L„^(x)] are of the same order (see 
(12.71) ). Since p > 2, by combining (I3.34p and (13.351) . we get 



^ _ var[L„^(x)lz3, 



var[L„^(x)] 



Since Ai(x, n^) < A(x, n^), the conclusion of Lemma [3.51 follows from 
(I3.28P by passing to the limit. □ 

The following result shows that the variance of effective noise and 
that of the its linear approximation are asymptotically equal 

Lemma 3.6. Let f and (^„) be as in Theorem \2.1\ Assume that the 
Lyapunov condition (12. 9p holds and that (ak) is a sequence of positive 
numbers that satisfy (I2.10p . Then, 

hm MM^,^k) -f-^{x))ls,] ^ ^ 

fc-oo cri var[L„^ (x)] 

where {Bj.} are the events defined in Section l3. 3. 1\ 
Proof. If (cTfc) satisfies (I2.10p . we know from Section [3.3. II that 
lim ||(w„,(x,cTfe) -/„,(x))lijj|2 = 

fe— >oo 

Therefore, by Lemma 13.51 it is enough to show that (13.280 holds for 
dk = ttk, where is defined by (I3.20p . 
Notice that 

^-2(p-i) {A{x,nk)f ^ ^ I II maxi<j<„.^, g(A(x, nfc))V fc 
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Thus, the conclusion of Lemma 13^ follows by passing to the limit. □ 

Theorem 3.7. Let f and {^n) be as in Theorem \2.1\ Assume that the 
Lyapunov condition ( \2.9\\ holds and that (cr^) is a sequence of positive 
numbers that satisfy (12.101) . Then, the process Wn^{x,ak) defined by 
fl2.13p converges in distribution to the standard Gaussian. 

Proof. We know from Theorem 13.31 that 0"^) converges in distri- 

bution to the standard Gaussian. Notice that 



^' var[(x„,(x,orfe) 

The conclusion of Theorem 13.71 follows from Lemma [3.61 and the con- 
verging together^'' Lemma. □ 

3.3.4. Noise with finite moments of order p > 4. Now, we will prove 
Theorem 12.11 under hypothesis H2. 

We will show that if the sequence has moments of order p > 4, 
then the choice of noise level (cxfc) (and hence the choice of outliers) is 
improved. 

Theorem 3.8. Let f be as in Theorem \2.1l assume that the Lyapunov 
condition (12. 9p holds. Let (^„) be a sequence of independent random 
variables with p > 4: finite moments, and that it satisfies A3(p). 

For any sequence (ak) that satisfies (12. lip , let {Bk} be the events 
defined by (I3.17P . Then, the process 

^ _ (xn(x,afc) -/"''(3;))lgJ 

a^ar[L„,(a:)] ^''''^ 

converges to the standard Gaussian. 

Proof. It suffices to verify that 

= ^iiZii 

satisfy condition (13.281) in Lemma [3.51 

Since ^y A2{x, Uk) < Ai{x, Uk) < A(x, Uk), we have that 

^^2(p-i) {A{x,nk)f f II maxi<j<„, \^j\\\l{A{x,nk)fak ^ ^^'^ 

A2(a:,nfc) " A2{x, Uk) 

The conclusion of Theorem 13.81 follows from (12.1 ID □ 

Remark 3.3. Notice from the definition of the outliers, that if the 
sequence of noises (^„) is supported on an interval [—b,b], then for k 
large enough Bk = fl, that is, the outliers are empty sets. 
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3.4. Proof of Berry— Esseen estimates (Theorem 12.21) . In this 
section, we present a proof of Theorem 12.21 First, we will prove two 
results (Lemmas 13.91 and 13. lOp that show how much the distribution of 
a random variable changes when we add a small perturbation. 

Lemma 3.9. Let ^ he a real random variable defined on some proba- 
bility space {Q,J^,¥), let n be a probability measure on (M, S(M)). Let 
h be a bounded Borel measurable function in M, and B & T. Then 

|E [h {ils)] - ^ih\ < 2\\h\\^¥[B^] + |E [h{i)] - i^h\ (3.38) 

where = J h{x)fi{dx). 

Proof. From the identity 

we have 

|E[MO]-E[/i(as)]| <2||/i||ooP[5l 
We obtain (13.381) by the triangle inequality □ 

Lemma 3.10. Let ^ and rj be real random variables defined in some 
probability space {fl, JF, P). Assume that for some Lipschitz function G 
on M and some constant < a < 1 we have 

1) snp,^^\F[7] < z] - G{z)\ < a 

2) E[\^-r^\]<a' 

Then 

sup |P[^ <z]- G{z)\ < (2 + L)a 
where L = L{G) is the Lipschitz constant of G. 

Proof. Define the event A = [\C, — r]\ < a]. Then by Chebyshev's 
inequality, we have ¥[A'^] < a. 

Partitioning the probability space Q, we have 

<z]= P[{^ <z}nA]+ P[{^ <z}n A'] 

< ¥[ri < z + a\+a 

Adding and subtracting G{z + a) — G{z) we obtain 

sup{F[^< z]~G{z)} < {2 + L) a (3.39) 

Similarly, by partitioning the space we have 

F[r]<z-a] = r[{r] < z - a} n A] + P[{r/ < z - a} n A""] 
< r[^<z]+a 
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Adding and subtracting G{z) — G{z — a) we obtain 

- a{2 + L) < M{¥[^ < z] - G{z)} (3.40) 

We finish the proof by combining (13.391) and (I3.40p □ 

Notice that from Holder's inequahty, if the Lyapunov condition (12.91) 
holds at a point x E M for some number 2 < p, then it also holds for 
s = min(3,p). By the same token, if (^„) is a sequence of independent 
random variables with p finite moments that satisfies A3(p), then it 
satisfies A3(s). 

If we denote by 

_ As{x,nk) 

and let bk be the left-hand side of (12.151) . The assumption (12.151) is 
expressed as 

bk < al 

Recall the processes Wn^. and In,, defined in (I2.12p and (13. 3p respec- 
tively. By Theorem 12. H we know that there are events such that 

Using Lemma [3.91 with h = l(_oo,z] and fi{dt) = ^'{t)dt, we get 
sup |P [Iklsk <z]- ^z)\ < 2bk + ak< fa^ 

for k large enough. Then, the conclusion of Theorem 12.21 follows from 
Lemma [3.101 by letting ^ = w^Xski ^ = luk^Bki Cl{z) = ^{z) and 
a = 3afc/2. □ 

Remark 3.4. Notice that the sequence of sizes of noises {(Jk}k and 
the sequence of outliers {Bk}k depend on the initial condition xq of the 
process Xn- Therefore, the central limit theorem in Theorem 12.11 and 
the Berry-Esseen estimates in Theorem 12.21 are point wise results. In 
section 13. 6[ we will give an example where the convergence to central 
limit theorem holds uniformly with respect all initial conditions. 

3.5. Cumulants. The Lyapunov functions Ag used in Theorem 12.11 
have a simple interpretation in terms of rather well known statistical 
quantities called cumulants in the statistical literature or Wick-ordered 
moments in statistical mechanics. 

The cumulants, for a given random variable ^, are formally defined 
as the coefficients Kp in the power series 

°° K 

\ogE[e''^]=J2-r('tr (3-41) 

p=0 
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Notice that i^sl^] = var[^]. It can be seen |Pet75l Ch. IV] that Kp[^] 
are homogeneous polynomials in moments, of degree p. Furthermore, 
for any a G M, and any pair of independent random variables ^, rj we 
have that 

Thecumulantsir3[e]/(^2[e])'/'and K,[^]/{K2[Q^ of^ = ^/./K^], 
called skewness and kurtosis respectively, are empirical measures of 
resemblance to Gaussian: the closer to zero the closer ^ is to a Gaussian 
|MGB74l p. 78]. Observe that these measures are scale invariant. 

Definition 3.1. For each positive integer p, we consider the cumulant 
functional Ap defined by 



n 



(3.43) 



Observe that |Ap(x,?T,)| < Ap{x,n) and that if p even, Ap = Ap. 
Furthermore, the cumulant functionals Ap{x, n) satisfy the analog to 
f l2.5l) for the Lyapunov functionals, namely 

Ap(x,n + m) = ((n'or(x))^Ap(x,n)+Ap(r(x),m) (3.44) 

In section 14.51 we will use (13.441) to study the asymptotic behavior of 
the cumulants in the case of systems near the accumulation of period 
doubling. 

Notice that when (^n)„ is an i.i.d sequence with p G N finite moments, 
we have that 



K,[Ln{x)] = Ks[^^]K{x,n) 



(3.45) 



for all < s < p. 

If ('Cn)n is a sequence of random variables that satisfies A3(p) for an 
integer p > 2, then from (12.41) and (I3.42p 



Ap{x,n) 



{A2{x,n)) 



p/2 



< Cr, 



Ap{x,n) 
{A2{x,n)) 



p/2 



(3.46) 



for some constant Cp > 0. Therefore, condition (12. 9p in Theorem 12.11 
can be interpreted as a measure of closeness to Gaussian. 

Using renormalization theory, we will obtain sharp asymptotics for 
(I3.46P for maps near the accumulation of period doubling and critical 
circle maps with golden mean rotation number. 
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3.6. Examples. In this section, we consider a nontrivial example of 
a system to which Theorem 12.11 applies. We also show that systems 
with enough hyperbolicity do not have a central limit theorem in the 
direction of Theorem 12. 1[ 

3.6.1. Smooth diffeomorphisms of the circle. We start this section by 
showing a central limit theorem for conjugate maps. We will assume 
that M = / or 

Lemma 3.11. Let f,g & C'^{M), and let (^„) he as in \2.1\ Assume 

f = h"^ o g o h 

for some map h such that h,h^^ G C^{M). The Lyapunov condition 
(12.91) holds for f at some point x if and only if it holds for g at h{x). 

Proof Notice that p = h'^ o g^ o h. Since h, h'^ e C^{M) and M 
compact, we have that 

c\ig^yohix)\ < \ {f^nx)\<C\ig^yohix)\ 

for some constants C,c> 0. If we denote y = f{x), we have 

c^Af(y,n)<A{(a;,n)<C^Af(y,n) 

for all < s < g. Now, Lemma [3.111 follows immediately. □ 

Recall that every orientation-preserving homeomorphism / of the 
circle is given by f{x) = F{x) mod 1, where F is a strictly increas- 
ing continuous function such that F{x) = f{x) + 1 for all x G M. For 
any a G M, we denote by Ra{x) = x + a mod 1 the rotation of the 
circle with rotation number a. 

We will use Lemma [3. Ill to show that Theorem 12.11 applies for smooth 
diffeomorphisms of the circle with Diophantine rotation number. This 
will be a consequence of the following Herman type theorem: 

Theorem 3.12. Let k > 2 and assume that f G C^(T^) is an orien- 
tation preserving diffeomorphism, whose rotation number a is a Dio- 
phantine number with exponent (3. That is 



for all I G Q, where C > is a fixed constant. If (3 < k — 2, then there 
is a map h G C''^^~^ such that 

f = o o h 
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Theorem 13.12^ due to |SK89j , is the most comprehensive global re- 
sult on the problem of smoothness of conjugacies of diffeomorphisms 
of the circle with rotations. The first global result on this problem 
is due to |Her79j for > 3, and generalized later by |Yoc84] . The 
first mayor improvement {k > 2) in the solution of this problem was 
obtained by p<S87l IKU891 [SK89] . A brief historical account of the 
problem and a detail proof of this and related results can be found 
in [SK89j . The technique developed in |SK89j constitutes one of the 
most important applications of the renormalization group method in 
dynamical systems. 

Example 1. Notice that Theorem 12.11 applied to is the classical 
central limit theorem. Thus, we can conclude that any Diophantine 
diffeomorphism of the circle has a Gaussian scaling limit for weak noises 
along the whole sequence of integers. Moreover, if the sequence of 
random variables S,n have finite moments of order p = 3 and we let 
cr„ = n"^/^, we can define the outliers by 



5^ 



II /"II Co max 1^/1 > Vn 



Notice that the outliers are independent of the initial point of the orbit. 
Furthermore, by Theorem 12.21 we have 

sup |PK(x) < ^] -$(2)1 < ^ 
(x,z)eTixR yn 

for some constant A> Q. □ 

3.6.2. Systems with positive Lyapunov exponents. In this section, we 
give an example of a system that does not satisfy the Lyapunov con- 
dition (12.91) . and for which the conclusion of Theorem 12. II fails. 

Example 2. Consider the map on R (or T-*^) given by 

f{x) = 2x 

Condition ( 12. 9p is not satisfied, and as we will see the conclusion of 
Theorem l2.1l fails. We will show that there is a limit for the scaled noise, 
not necessarily Gaussian, which depends strongly on the distribution 
of the original noise Notice that 



Wr. 



^nvar (e;^, 

2V1 -4-" ^ ^ 
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If is an i.i.d. sequence with uniform distribution U[—l, 1] then, Wn 
will converge in law to a compactly supported random variable with 
characteristic function 



n 

k=2 



2y2sin(2"'=V23^) 



On the other hand, if {^n) is an i.i.d sequence with standard normal 
distribution, then Wn has standard normal distribution for all times. 

Similar results hold for hyperbolic orbits and the reason is that 
derivatives at hyperbolic points grow (or decay) exponentially, that 



P 



for some number a > 0. If is an i.i.d 



sequence for instance, then the cumulants Kp of order p > 2 oi the 
normalized variable L„ do not decay to zero. In fact we have 



LJx) 



A/var[L„(x)] 



where Kp is the cumulant of order p of ^i. Therefore, for hyperbolic 
orbits the scaling limit depends strongly on the distribution of the 
sequence □ 
In view of the example above, we can see that the scaling version 
of the central limit theorem we have proved does not apply for maps 
which are hyperbol ic. In that res pect, it is curious to mention that 
a celebrated result GS96| |Lyu97[ IGS97 shows that the maps in the 



quadratic family are hyperbolic for a dense set of parameters. 



4. Central limit theorem for maps the accumulation of 

period doubling 

In this section, we will show that the orbits (11.11) at the accumula- 
tion of period-doubling satisfy a central limit theorem (Theorem 12. 3p . 
We consider some statistical characteristics of the effective noise, i.e. 
Lindeberg-Lyapunov sums, and show that they satisfy some scaling 
relations (Corollary 14. 8p . We use this scaling relations to show that 
the Lyapunov condition 12.91 holds for certain class of initial conditions 
(Theorem I4.15P . 

4.1. The period doubling renormalization group operator. We 

will consider analytic even maps / of the interval / = [—1,1] into itself 
such that 

PI. /(0) = 1 

P2. xf'ix) < for X ^ 
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P3. Sf{x) < for all a; 7^ where 5* is the Schwartzian derivative 
{Sf{x) = {fix)/ fix)) - (3/2)(r(a:)//'(x))2) 

The set of functions that satisfy conditions P1-P3 is called the set of 
unimodal maps. 

Let us denote by A/ = /(I) and 6/ = /(A/). We will further assume 
that 

P4. < |A/| < hf 
P5. < f{hf) < \Xf\ 

The set V of functions that satisfy P1-P5 is called the space of period- 
doubling renormalizable functions. 

It follows from P4 that the intervals Jq = [— |A/|, |A/|] and Ji = [&/, 1] 
do not overlap. If / G P, then /o f\i^ has a single critical point, which 
is a minimum. The change of variables x ^ XfX replaces Jq by I and 
the minimum by a maximum. 

Definition 4.1. The period doubling renormalization operator T on 
V is defined by 

Tf{x) = ^fof{\fx) 

If f,Tf,..., T^^'-^^f e V, we denote by r„ = f"{0) (or if we need 
to emphasize the dependence on /). We have that 



T-fix) 



A 



Tnf 



n+1 



The renormalization group operator T has been well studied by many 
authors, jMTTl [TC781 ICEgOl ICELSOl ILan82l IVSKM IEW851 [E^s86 



[CL891 [Sul92l IMar98l IdMvSQSl IJS02j . The following results, whose 



proofs can be found in the reference above, will be very useful for our 
purposes. 

Rl. For each e N, there is an analytic function hk such that 
gk{x) = hk{x'^^) is analytic on some open bounded domain Dk C 
C containing the interval [—1, 1], gk restricted to I is concave, 
and for all x G / 

Tgk{,x) = gkix) (4.1) 

Furthermore, since Sf{x) < for all a: 7^ 0, it follows that 
g"{x) < for all x 0. 

Notice that (14.11) implies that 

g'ki^) = g'kigk{hx))g'kiXkx) 
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Then, by taking the hmit a; — we get 



R2. The domain Dk can be chosen so that dD^ is smooth, C 
Dk and gtih Dk) C Dk. 
This condition is closely related to the fact that is in the domain of 
the period doubling renormalization operator. 

For each k, denote by H{Dk) the space of analytic functions on D^, 
and let Qk C H{Dk) be the Banach space of analytic functions 

Gk = {f: f{z) = z^'fiz), f{z) e HiDk), f e Com} 

endowed with the sup norm. Notice that if f e Gk, then f^^\0) = 
for all derivatives of order < j < 2k. 

R3. There is a neighborhood V of g^ in {Q^. + 1) fl P where T is 
differentiable. 

R4. A consequence of R2 is that if / G V, the derivative DT{f) is 

a compact operator on . 
R5. T admits an invariant stable manifold Ws((7fc); such that if / G 

Ws{gk) then, 

hm \\T-f-gk\\ciD,) = 0{ul) (4.2) 

where < cj^ < 1 is a universal constant. 

The set Ws{gk) is called the universality class of gu. 

If we denote by A^,^ = A^ then, for each G N we have that V^^ = 
(Afc)"'. Moreover, the exponential convergence (R4) implies that for any 
/ G yVs{gk), At" converges to Afc exponentially fast. 

4.2. Renormalization theory of the noise. In this section, we de- 
velop a theory of renormalization for the noise. Since the renormaliza- 
tion group operator T gives information for times 2" at small scales, 
we will study first the effective noise at times 2" for orbits x„ starting 
at zero. 

Let / G Ws((?) be fixed, and for each p > denote by k^\x) = 
A^ix, 2"). From we have 

"'A;("-i)(a;) + 4"-i)o/2""'(x) (4.3) 



kir\x) 



Since f E V, then —f\x) > 1 > for all a; G [/(A/), 1]. Therefore, by 
taking a determination of logarithm defined on 

C \ {x + : X < 0, y = 0} 
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we can analytically define the map z ^ (— /' o /(-^z-^))^ on -D/;. In 
particular, notice that if 2; G -0^ fl M, then 

(-/o/(A;^))^ = |f o/(A^^)r (4.4) 

4.2.1. Lindeberg-Lyapunov operators. The following family of opera- 
tors will be very useful in the study of the propagation of noise. 

Definition 4.2. For each / G Wsig), the family of Lindeberg-Lyapunov 
operators {/C/,p '■ P > 0} acting on the space H'''{Dk) of real and 
bounded analytic functions on D/. which are continuous in D^, is de- 
fined by 

^/.^'^W = ° h{Xfz) + hifiXfz))] (4.5) 

where h e H\{Dk). 

The operators defined by (14. 5 p belong to a class of operators called 
transfer operators, which have been extensively studied [BalOOj May80| . 
It follows from (gSD and (gS]) that if 

X = T^z, = |r{| k^p^iVl^z), 

then, defining the function 1 by 1(2;) = 1 

fcJ")(z) = (/C^„-.;,,---%pl)(^) (4.6) 

Equation (14.61) provides a relation between the Lindeberg-Lyapunov 
operators defined by (14. 5 p and the effective noise. Indeed, recall that 
the sequence of noises (^n)n satisfies < inf.„ ||^||2 < sup„ ||^||p < c. 
Then, for s = 2, p we have that 

2" 

This means that (14. 6 p estimates the growth of the linearized propaga- 
tion of noise at times 2" with initial condition at the origin. 

4.2.2. Exponential convergence of the Lindeberg-Lyapunov operators. 
Recall that for a map / e Wsidk), we have that T^f{z) converges 
exponentially fast to gk- Our next result implies that the sequence 
of Lindeberg-Lyapunov operators /Ct" /,p also converges exponentially 
to )Cgf,,p- This implication will be important in Section 14. 3[ Proposi- 
tion 14.71 where we analyze the asymptotic behavior of the product of 
operators in (14. 6p . 
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Lemma 4.1. Let f G yVs{gk) close enough to Qk- For all p > 0, /C/,p 
is a compact operator on H'^{Dj^). 

Furthermore, there is a constant Cp such that 

\\^f,P - ^9k,p\\ < ^pWf - 9k\\c{Dk) (4-7) 
Proof. For any set D C C, let = {z : d{x, D) < e}. Recall that 

AfcA. C Dk 
9k{hF>k) C Dk 

Since for all z E 

\\fZ-Xkz\ < \\f - gk\\r),,diam{Dk), 

we can choose a neighborhood U of gk such that for all / e VVsigk) 

XjWk C C Dk (4.8) 

where e > is small enough. Let us denote by Vk — XkDk . 
Notice that for all z G -Dfc, 

\gk{>^kz) - gk{\fz)\ < diam(i:>fe)||/-^fe||i3, 

\f{\fz)-gk{Xfz)\ < \\f-gk\\vk 

Then, for all z E Dk we have that 

\f{Xfz)-gk{Xkz)\<C\\f-gk\\D, (4.9) 

for some constant C. Therefore, shrinking U if necessary, we can as- 
sume that 

fi^fDk) C gkiXkDk)' C Dk (4.10) 

Let us denote by Wk = gk{>^kF>kT- 
For all z e Dk we have that 

\f'{f{Xfz))-g'k{gkiXkz))\\ < \g'kigkiXk z)) - g'k{f (Xk z))\ 

+ \g'k{f{Xfz))-f{f{Xfz))\ 
Using Cauchy estimates and the following inequalities 

Ig'kigkihz)) - gk{f{Xkz))\ < Wg'kWwkWgk - fWvk 
l^?U/(A/^))-/'(/(A/;^))| < \\g'k-f\\w,, 
we obtain that 

II/' ° f°Xf-gk° 9k o Afcll^^ < C\\gk - /Ibfe (4.11) 

Let p > be fixed. Notice that the Lindeberg-Lyapunov operator /C f^p 
has the form 

lCf,ph — Ufh + Rfh 
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where U : h \—>- [tp o f o f o Xj)[h o Xj) and R\ h^hofoXj^ where 
il){z) = is defined onO = C\{x + z?/:a;<0, ?/ = 0}. By shrinking 
U if necessary, we can assume by (14.111) that if / G Wsigk) then, 



Let us denote by Yk = g'k{9k{^kDk)) From fl4.9p . we have that 

(4.12) 

< C\\h\\Djf - gkWo, 

Prom (14. lip , we have that 

||(t//-f/,J/^||<C'||^||yJ|/.||5j|/-^7.||s, 

+ C\\^'\\Yjh\\vJf-gk\\5, (4.13) 
The exponential convergence follows from (14.120 and (14.130 . 

Notice that the operators JCf^p send bounded sets in H^{Dk) into 
bounded sets in W{Dk). 

Furthermore, for all Zi, Z2 E Dk we have from (14. 8 p and (14.101) that 

\Rfhizl) - Rfh{z2)\ < \\ihofy\\yJzi-Z2\ 
\Ufh{z,)-Ufh{z2)\ < \i{ijof)hy\^Jz,-Z2\ 

The compactness of the operators ICf^p follows from Cauchy estimates 
and the Arzela-Ascoli theorem. □ 

We will study the spectral properties of the Lindeberg-Lyapunov 
operators to show that (12. 9p holds, and then by Theorem 12.11 we will 
obtain the Gaussian scaling limit for systems at the accumulation of 
period-doubling. 

4.3. Spectral analysis of the Lindeberg-Lyapunov operators. 

Let C be the cone of functions on H^{Dk) which are non-negative when 
restricted to Dj. fl M. Note that intC, the interior of C consists of 
the functions in H^{Dk) which are strictly possible when restricted to 

Let / G yVs{g) and p > fixed, and denote JCp = ICj^p. 
We have that 

(i) x:,(c\{o})cc\{o} 

(ii) ;Cp(int(C)) C int(C), and 

(iii) for each /i G C \ {0}, there is an integer n = n{h) such that 
la^h G int(C). 
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Properties (i), (ii) are quite obvious. Property (iii) follows from the 
observation that if a function is non-negative and strictly positive in 
an interval /* - the set of zeros has to be isolated - then, by the 
definition of K,p (14 .Sp it is strictly positive in the image of /* under 
the inverse of the maps A/-, /(A/-)). Since the functions A/, /(A/-) are 
have derivatives strictly smaller than 1, it follows that a finite number 
of iterations of the inverse functions covers the whole interval. 

By Krein-Rutman's Theorem |Sch71t p. 265], |Tak94j . we have that 

Proposition 4.2. Let f G Wsig) and p > 0, and denote by JCp = /C/,p- 
Then, 

Kl. The spectral radius pp of JCp is a positive simple eigenvalue of 
ICp. 

K2. An eigenvector ipp G W'^Dk) \ {0} associated with pp can he 
chosen in int{C). 

K3. // fi is in the spectrum of ICp, fi pp, then fi is an eigen- 
value of ICf^p and < Pp. 

K4. If h & C \ {0} is an eigenvector of JCp, then the corresponding 
eigenvalue is Pp. 

Remark 4.1. We emphasize that Krein-Rutman's theorem only needs 
to assume that the operators are compact and that preserve a cone. 
Hence, the only properties of the operators Kp that are used are pre- 
cisely i),ii),iii) above. Later in Section Owe will see another applica- 
tion of the abstract set up. In the case of Section O the space will be 
a space whose elements are pairs of functions. It is important to note 
that Corollary 14.31 Proposition 14.41 and Proposition 14.61 remain valid in 
the case of several components. 

We will say that two sequences of functions {a„(z)}„ and {bn{z)}n 
of function in H^{Dk) are asymptotically similar (denoted by a„ x 6„) 
if 

lim Wan/bnWurfD^) = 1 

As standard, we use the notation ip > (p to denote ip — (f) & intC. In 
particular > means G intC. For function in H^{Dk), (p > ip is 
just means that (p{z) > ip{z) for all z E DkH M. 

One property of the cones in the spaces we are considering in this 
Section and in Section [5] is that: 

property 4.1. Let (p,ip E intC Then, there exists 6 > such that 
^p > 5<p. 



A consequence of Proposition 14.21 is: 
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Corollary 4.3. h E C \ {0} then, there is a constant Cp = Cp{h) > 
such that 

\\}C;h{z) - Cpp;M^)\\ < C{pp - 6r (4.14) 

for some C > 0, 5 > 0. 

In the one- dimensional case, this implies. 

lC;h{z) X c,p;V^,(z) (4.15) 

To prove ( 14.14p we note that we can define spectral projections corre- 
sponding to the spectral radius and to the complement of the spectrum. 
\{ h = Cpipp + h"^ with h"^ in the space corresponding to the spectrum 
away from the spectral radius. 

By the spectral radius formula, we have ||/Cp/i^|| < C{pp — 5)" for 
some 6 > 0. On the other hand, we have that for some S > 0, h > 5ipp 
over the reals. Hence, IC^h > Sppipp over the reals. By comparing with 
JCph = Spptjjp + /Cp/i<, we obtain that Cp > 0. 

The conclusion (14.151) follows from the fact that ipp is bounded away 
from zero. 

Another important consequence of Proposition 14.21 and the analyt- 
icity improving is that the spectrum and the eigenvalues are largely 
independent of the domain of the functions we are considering. This 
justifies that we have not included it in the notation. 

Proposition 4.4. Assume that Dk, Dk are domains satisfying the as- 
sumptions i?2- 

Assume that applying a finite number of times the inverse functions 
of either Xk- org{Xk-) we obtain a domain which contains D^. Then, 
the spectrum of JCf^p acting on H'^{Dk) is contained in the spectrum of 
JCf^p acting on H^{Dk). 

We know that the spectrum is just eigenvalues (and zero). If we use 
the functional equation satisfied by an eigenvalue of /C/^p, we see that 
the eigenfunction extends to the images of D^. under the inverse images 
of Afc- or g{Xk-) 

Hence, by the assumption, we obtain that the eigenvalues in H{Dk) 
are eigenvalues in H{Dk). □ 

Of course, if the hypothesis of Proposition 14.41 is satisfied also when 
we exchange Dk and Dk, we conclude that the spectrum is the same. 

Note that in our case, since the inverse maps are expansive. Proposi- 
tion 14.41 shows that the spectrum of the operator is independent of the 
domain, provided that the domain is not too far away from the unit 
interval. 
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For the case of the quadratic fixed point there is a very detailed study 
of the maximal domains in [ELSlj . The results of [EL81j imply that 
there is a natural boundary for the eigenfunctions of the operators. Of 
course, if we consider domains larger than that, the results are different. 
Nevertheless, for domains inside this natural boundary, the spectrum 
is independent of the domain. 

4.3.1. Properties of the spectral radius of the Lindeberg-Lyapunov op- 
erators. For / G yVs{g), we have that the spectral radii of the operators 
JCf^p satisfy some convexity properties. 

Theorem 4.5. Let f be a map in Ws{gk) close enough to g^. For each 
p>0, denote by ICp = ICf^p, and the spectral radius of Kp by Pp. Then 

SI. 

(A7V'(l)r <Pr>< (A7V'(l)r + (-A/)-^ (4.16) 
In particular, if f = gu then 

1< Af ^Pp < 1 + (-A,)(2'=-i)P (4.17) 

52. For p > 0. The map p ^ pp is strictly increasing and strictly 
log-convex. 

53. The map p \—>- \og[pp\/p is strictly decreasing. 

For future purposes in this paper, [S3.] is the most important of the 
consequences claimed in Theorem 14.51 

Proof. SI: For fixed p > 0, let ipp be a positive eigenfunction of JCp. 
Denote by 6 = (— A/)"^, and notice that b > 1. Let Up, R and Sp be 
the nonnegative compact operators defined by 

Uph{z) = {-fof{Xfz)rh{Xfz) 
Rh{z) = hof{\fz) 
Sph{z) = {-f{l)Yh{Xfz) 

Since K,ph = h {Up + R) h, for any /i G C \ {0} and any n G N we have 
that, by the binomial theorem 

/C;/i > 6"(?7p" + i?")/i > 6" (4.18) 

Notice that 

If we take h{z) = ipp{z) and then let z = 0, the left-hand side of (14.161) 
follows. 
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Since 1 < ~f'{z) < -/'(I) for all z G [/(A/), 1], we have that 

}C;h < V\Sp + Rfh 

positive function h. Notice that the spectral radius of the positive 
operator S'p + i? is (— /'(l))^ + 1. Then, the right-hand side of fl4.16l) 
follows. 

The claim in [SI.] for the special case f = Qk follows from 

To prove [S2.], [S3.] we will use the following Proposition 

Proposition 4.6. In the conditions of Proposition Assume that 
(f) is a positive function X is a positive number such that 

K,p(j) > A0 

Then, 

Pp> X 

Note that because the functions are strictly positive, we have /Cp0 > 
(A + for some 6 > 0. Then, applying this repeatedly, we obtain 
}C^(f) > (A+5)"(^. By CoroUaryiJl we have }C^(j) = Cpp;jVp+o((p-5)"). 
Since the function ipp is strictly positive, we have the desired result. □ 

Now, we continue with the proof of [S2.], [S3.] 

We recall that the operator JCp has the structure: 

}Cph = A^ho f\ + ho f2 (4.19) 

where A{z) = —f'of{Xfz), fi{z) = XfZ, f2{z) = f{Xfz). In particular, 
A>1. 

The strict monotonicity of Pp follows from the fact that if g > p 

Pqi^q = A'^lpg O /i + O /2 
> APtfjg O /i + 7/;^ O /2 

To prove that log(pp)/p is strictly decreasing, we will show that for 
any < p and any a > 1 we have ppa < Pp, which is equivalent to the 
desired result. 

If we raise the eigenvalue equation for }Cp to the a power, we obtain, 
p^^; = iAP^Pp o A + ^p o f,r 

> + i^pOf2 (4.20) 

The inequality above is a consequence of the binomial theorem for frac- 
tional powers. Applying Proposition 14.61 we obtain the desired result 

Ppa < Pp ■ 
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To prove the strict log-convexity of pp, we argue similarly. We will 

1/2 1/2 

show that for < p < g and for p(p+g)/2 < Pp Pq , which implies strict 
convexity. 

Multiplying the eigenvalue equations satisfied by the spectral radius 
for /Cp and for /Cg and raising them to the 1/2, we have: 

> (AP+''(^p^,)o/i + (^p^,)o/i)^/' 

> A(^+^)/2(^p^,)V2 o A + (^p^,)V2 o 

The first inequality comes by expanding the product and ignoring the 
cross terms, and the second inequahty is a consequence of (1 + x)^^"^ < 
1 + x^/^ for a; > 0. 
Applying Proposition 14. 6[ we obtain the desired result. □ 

□ 

4.4. Asymptotic properties of the renormalization. From Lemma 
fl4.ll) . we know that \\T'"'f — gkWok = ^(^fe) implies 

II^T"/,p-/C,„p||=OK) (4.21) 

for all p > 0. Therefore, the asymptotic properties of JCT"f,p become 
similar to those of JCg^,p for large n. In Proposition 14.71 below, we 
show that the dominant eigenspaces of the Lindeberg-Lyapunov oper- 
ators, after a large number of renormalizations, align to the dominant 
eigenspace of the cumulant operator at the Feigenbaum fixed point gk- 

Proposition 4.7. Let pp be the spectral radius of /Cp = /Cgj.,p. There 
exists a neighborhood U of gk such that if f E Ws{gk) and h is a 
positive function h G then, 

c-'p; < ICrn-i /,p ■ ■ ■ lCf,ph{z) < cp; (4.22) 

for some constant c. 

The proof of this Proposition will be done after stating a conse- 
quences. A Corollary of (14.221) is the following 

Corollary 4.8. For any p > 0, let Pp be the spectral radius of /Cp = 
1. There are constants < c < C such that 

A^(rj;z,2") A^(rj;z,2") 

c < inf " ' / < sup r I r 1 < C 4.23 
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2. If f E yVs{gk) is dose to Qk then, for any p > 2 

A,^(r{z,2") 

lim sup — J — J = U 

z&i {A^(r^2,2")}p/2 

Proof. (1.) From (14 .Op and Proposition 14.71 we have that 

cp;\ri\p<A,{ri^X)<Cp;\rir 

for constants < c < C. Part (1) follows from the exponential conver- 
gence since (14 ■2p implies that 

c\\kr<\ri\<c\\kr 

for some positive constants c and C. 

Parts (2) follows from part (1) and Proposition 14.51 since for all p > 2, 
we have pp < a/pI ■ ^ 

The proof of Proposition 14.71 is based on the following general result 
(Proposition (14. 9p ) on convergent sequences of operators on a Banach 
space. The method that we will use is a standard technique in hy- 
perbolic theory |HP70j . where it is used to establish the stability of 
hyperbolic splittings 

Lemma 4.9. Let X be a Banach space, {/Cn}^i be a convergent se- 
quence of operators on X such that lim„ ||/C — /C„|| = 0. Assume that 
there is a decomposition of X = © , and constants < A_ < A+ 
such that E"^'^ are closed subspaces and 

IC{E>) = E> 

1C{E<) C E< 

Spec(/C|B<) C {^gC:|^|<A_} 

Spec(/C~^|s>) C |2GC:|z|<^|. 

There exist e > such that if sup„ ||/C,i — /C|| < e, then there are 
sequences of linear operators {A„}5^q C C{E'^,E^) and {_B„}^q C 
C{E^, E^) with livcin \\An\\ = and lim„ ||-Bn|| = 0, such that the spaces 

E> = {x®Ar,x: xeE<} (4.24) 

E< = {x®BnX: xeE<} (4.25) 

satisfy X = E^ © E^ and 

E> = K,n---K,i{E>) (4.26) 

E< = /C„---/Ci(Eo<) (4.27) 
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Proof. Since X is the direct sum of invariant closed subspaces with 
respect to /C, using adapted norms |CFdlL03| Appen. A] we can assume 
without loss of generality that 

\\]C\e<\\ < A_ (4.28) 

\\lC-%>\\ < ^ (4.29) 
^+ 

and that if x = + with x*^ G and x^ G , then ||x|| = 
||x^|| + llx-^ll. 

The operators /C and JCn are represented by the following matrices 
of operators with respect to the fixed decomposition X = E^ © E^ 

K = (4.30) 

K„ = (^,t> Jl,>>) (4.31) 
where ||i;<|| + \\S<>\\ + ||i><|| + - 0. 

I. We show first the existence of the sequence of transformations {An} 
in (Kmf . 

If use f l4.30p to compute the image of a point x + AnX we obtain that 
it has its E^jE"^ components are, respectively: 

{K> + 6»)AnX + 6<>x (4.32) 
{K< + 6«)x + An6> < X (4.33) 

The invariance condition that Kn{x + A„x) is in the graph of An^i 
is equivalent to saying that applying An+i to the second row of fl4.32p . 
we obtain the first. 

Therefore, the invariance condition is equivalent to 

A^+,[iK< + S«) + AJ> <] 

= {K> + (5»)A„ + 5<> 

Isolating the first An in the RHS, we obtain that the invariance 
fl4.26p for the spaces is equivalent to the following invariance equation 
for {An} 

An = (/C> + 5>>y^ {An^i (/C< + 5<<) + An+i5<>An - 5><} (4.34) 

Notice that if {A„} G A satisfies fOIjl . then i^[?M holds. 

We will show that there is a unique solution {An} of the invariance 
equation fl4.34p in an appropriate space of sequences. To be more 
precise, let A be the space of sequences of operators in C{E^, E^) that 
converge to zero. Observe that ||y4|| = sup„ ||A„|| defines a norm on A, 
which makes it into a Banach space. 
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Let T be the operator on A defined by the right-hand side of equation 
( 14341) . For A and A in ^ we have that 

\t{AU - r{AU <m> + S>>)-\A^+, - I„+i)(/C< + S<<)\ + 

|(/C> + 5>>)~'I|A„+A<>(A.-I„)| 

Let R > max(||yl||, \\A\\). Since 6n 0, we can assume with no loss of 
generality that for a given e > 0, sup„ ||5„|| < e. 

\\t{A) - r(I)|| < (A;U_ + t)\\A - A\\ + 2(A;^ + t)t\\A - A\\R 

Since A^1^A_ < 1, we can choose e > small enough so that 

e + 2i?(A;ie)e < 
(A;iA_ + e)/2+(A;^e)(e/?2 + e) < R 

we get that r is a contraction on the ball B{0;R) C A. Therefore, 
there exists a unique solution, A, to the invariance equation (I4.34p . 

IL An argument very similar to the one above shows that (14.271) holds if 
and only if {-B„} C C{E^ ,E'^) satisfies the following invariance equa- 
tion 

= {(/C< + (5<<)5„ + 5<< - 5„+i5><fi4 (/C> + (5>>)-i (4.35) 

Let B be the space of sequences in C{E^,E'^) that converge to zero 
with norm = sup„ \\Bn\\, and let f be the operator on B defined by 
the right hand side of (I4.35p . An idea similar to the one used in part (I), 
shows that f is a contraction on an appropriate closed ball. Therefore, 
there is a unique solution to the invariance equation (14.351) . □ 

Remark 4.2. \\An\\ and ||-Bn|| are measures of distance between the 
spaces E^'^ and E^'^ . Lemma IT9] states that these respective spaces 
align in the limit. 

The next Corollary gives control on the growth of the products of 
operators /C„. 

Corollary 4.10. Under the hypothesis of Proposition let us as- 
sume that dim{E^) = 1. Then: 

(1) For any v & X 

n n 

WiK - e,) \\v\\ < ||/C„ ■ ■ ■ ICovW < J](A+ + e,) ||^;|| (4.36) 

i=i i=i 

where e„ = e„(||/Cj — /C||), is such that Imit^Qenit) = 0. 
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(2) // ||/Cn — /C|| = 0(c<j") with < c<j < 1 , then we have that 

c-^Xl\\v\\ < \\ICn---}Cov\\ <cXl\\v\\ (4.37) 
for some constant c > 0. 

Proof. Using adapted norms, we can assume that fl4.28l) . fl4.29p hold, 
and that if x = x^+a;^ with x"^ e E'^ and G , then = ||a;^ || + 
\\x>\\. For s ="<", " >", let E^ = E' and define E^^^ = JCn{E'J. 
(1) Let V G E^, then by Lemma 1^9) there exists a unique w G E^ such 
that 

V = W + AnW 

Since ||v4„||, \\Bn\\ 0, we may further assume that 

sup II A„ II < 



sup||5„|| < 



2||/C|| 
5 



ll^n-/C|| < ^ 

for (5 > small. Since ||/Cw|| < A_||tf || < A+HtiH, from triangle inequal- 
ity we have 

< -\\V\\ + + Ar,w)\\ 

< iK + S)\\v\\. 

Let V G E^ then, by Lemma [4.91 there is unique z G E^ such that 

V = BnZ + Z. 

Therefore, ||f || < gl^lkll + Ikll' which is equivalent to 

'\v\\ 



\z\\ > 



1 + 6 

From the triangle inequality, we have that 

||A^n^^|| > ll^"^!! ~ II (^n. ~ 

> ||/C(5„z + ^)||-^||t;|| 

^ lir II ^l II 
II II 2" II 

> {K-S)\\v\\ 

Combining these results, we obtain inequality (14.36^ . 
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(2) In the case of exponential convergence, (I4.37P follows from an 
elementary result on sums and products, namely: if lim„ sup^- |c„j| = 0, 
and sup„ Yl]=i l^njl < oo, then 

+ Cnj) ^ if ^ Cnj c 

j j 
See IDurOSl p. 78]. □ 



Proof of Proposition |^. 7[ By Proposition 14.21 we have that that the 
spectral radius pp of Kp is a simple eigenvalue. Then, (14.221) follows 
from Corollary 14 . 1 01 and the exponential convergence, see (I4.2ip . of the 
Lindeberg-Lyapunov operators. □ 

Remark 4.3. The estimates (I4.23p for p = 2 are obtained in |VSK84] 
using the Thermodynamic formalism. 

4.5. Cumulant operators. In Section [375] we introduced the concept 
of cumulants, see (I3.4ip . One important fact is the cumulants of ran- 
dom variables, normalized to have variance one, provide an empirical 
measures of resemblance to Gaussian. In this section, we make a con- 
nection between the spectral properties of the Lindeberg-Lyapunov 
operators studied in Section 14.31 and the cumulants of the propagation 
of noise. 

Recall that the functions / G Wsigk) are defined in a domain Dk, 
see R2 in section HTTl 

Definition 4.3. For each / G yVs{g), the family of cumulant operators 
{/C/,p : p G N} acting on the space H^{Dk) are defined by 

= j^p [if ° /(A/^))" h{Xfz) + hifiXfz))] (4.38) 

By an argument identical to that given in Lemma W7L\ it follows from 
the exponential convergence of T"/ that 

Lemma 4.11. Let f G Wsigk) close to Qk- For each nonnegative 
integer p, JCf^p is a compact operator on H^{Di^). Furthermore, and 
/Ct"-i/,p converges to ICgj^^p exponentially fast. 

Notice from Definition 14.31 that 

ICfj, = Kfr, p even 
^ ll^/,pll P odd 
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Denote the spectral radius of /C/,p and /C/,p by pf^p and p/,p respectively. 
Then, from (14.391) we have that pf^p = pf^p for all p even, and for all 
integers p > 2 

~Pf,p < Pf,p < {Pf,2f' (4.40) 
Let / G yVs{gk) be fixed, and recall from Section [375] the cumulant 
functional A. For x G /, let 

kn,p{x) = r„^Ap(r„a;, 2 ) 

Then, from (13.44p . we have that 

fcn,p(s)(x) = /CTn-i/,p- ■ ■/C/,pl(a;) (4.41) 

We have that from Corollary 14. 101 that for very h G H^{Dk) there is 
a constant c such that 

c^^Pp < \\ICTr^-if,p---ICf,ph\\ < cp^ 

Recalling the notation Kp[C,] for the p-th order cumulant of a random 
variable ^, we have from (I4.40p that 



Kp 



{Tl x) 



Jv8iT[L2n{rlx)]} \P2 



Pp \ KA^) , Q ^4 42) 



as n oo, where hg^^p is a positive eigenfunction of /Cg^.,^. From the 
compactness of the cumulant and Lindeberg-Lyapunov operators, the 
asymptotic expansions (I4.42p are improved by adding more terms that 
involve eigenfunctions of smaller eigenvalues |FdlL06j . That is, we start 
from the standard asymptotic expansion of powers of a linear compact 
operator: 

N 

JCl^pUiz) = p;hg,,pc{u) + J2p1c,{u)^,{z) + 0{pl^,) (4.43) 

j=i 

where pj are the Jordan blocks associated to the eigenvalues pj (in 
decreasing order of size), Cj{u) are the projections on this space and 
j{z) is a basis for the eigenspace. 

For the nonlinear renormalization operator, the theory of |FdlL06] 
shows that there are very similar expansions for the asymptotic noise. 
The expansion involves not only powers of the eigenvalues of K, but 
also powers of products of eigenvalues. 

These asymptotic approximations are very reminiscent to the higher 
order asymptotic expansions in the central limit theorem, such as Edge- 
worth expansions |Fel71l p. 535]. (The asymptotic expansions (14.431) 
corresponds to Edgeworth expansions of 2" terms). It seems possible 
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that one could also prove Edgeworth expansions for the asymptotics 
of the scaling limits but it looks like the powers appearing in the ex- 
pansions will be very different form the standard semi-integer powers 
and will be related to the spectrum of some of the operators giving the 
renormalization of the cumulants. 

4.5.1. Numerical conjectures on spectral properties of the cumulant op- 
erators and the Lindeberg-Lyapunov operators. If p is even, we know 
from Proposition 14. 21 that pj^p is a that single eigenvalue of Kf,p = ^f,p- 
If p is odd, we have by compactness that /C/,p has an eigenvalue 

/i/,p e{z eC:\z\= pf^p}. 

In the quadratic case {k = 1), numerical computations [DEdlLOB] 
suggest that 

CI. If Spec(/Cg,,p) \{0} = {pn,p}n with \pn,p\ < IfJ'n-iJ, then Pn,p ~ 
A" (Here, u„ ~ f„ means that lim„M„/f„ = 1). 

C2. Spec(/Cgj,p) \ {0} is asymptotically real, that is: /i„ € M for all 
n large enough. 

For all p odd 

C3. Pg^^p is an eigenvalue of /Cgi,p. The rest of Spec(/Cg^^p) is con- 
tained in the interior of the ball of radius Pgi,p. 

C4. A^^^ ~ Pgi,p ~ Pgi,p as p ^ CX). 

C5. Spec(/Cgi,p) \ {0} = {pn,p}, where \jln,p\ < Wn-i,p\ is asymptoti- 
cally real. 

C6. Pn,p ~ Pn,p ~ A". 

Conjectures C1-C6 are related to conjectures on behavior of the Perron- 
Frobenius operators of real analytic expanding maps, see [BalOOl MaySO 



Rug94| . The case p = 1 corresponds to the problem of reality of the 



spectrum of the linearized period doubling operator. It is observed in 
|CCR90] that the spectrum in this special case appears to be real, and 
that all the eigenvalues behave as A". 

4.6. Proof of the central limit theorem for systems in a domain 
of universality (Theorem 12.31) . 

4.6.1. Introduction. In this section, we prove the central limit theorem 
(Theorem 12.31) for maps / in the domain of universality Wsigk)- 

The first step is to prove a central limit theorem and Berry-Esseen 
estimates for ) with initial condition x in the orbit of 0. 

We will verify the Lyapunov condition (12.91) along the whole sequence 
of positive integers. Then, we apply the result of Theorems 12.11 and 12.21 

The main observation is that the renormalization theory developed 
in Section U]2] gives control on As{x, 2") for x sufficiently small (roughly 
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< |r„+i|), see Corollary I4.8[ Then, the main tool to obtain control 
of the effect of the noise on a segment of orbit, will be to decompose 
the segment into pieces whose lengths are powers of 2 and chosen so 
that the renormalization theory applies. 

This decomposition will be accomplished in Section I4.6.2I The main 
conclusion of this section will be that this decomposition is possible 
and that, the final effect over the whole interval is the sum of the 
binary blocks - which are understandable by renormalization affected 
by weights. See (14.451) . The weights measure how the effect of the 
dyadic block propagates to the end of the interval. 

The effect of the dyadic blocks in the final result depends on some- 
what delicate arguments. This follows from the the following heuristic 
argument. 

Note that when an orbit passes though zero, the derivatives become 
very small (hence the effect of the noise up to this time is erased). 
Hence, when an orbit includes extremely close passages through zero, 
the next iterations tend to behave as if they were short orbits (and, 
hence the noise is far from Gaussian). On the other hand, the close 
passage to zero are unavoidable (and they are the basis of the renor- 
malization). 

Hence, our final result will be some quantitative estimates that mea- 
sure the effect of the returns. The main ingredients is how close the 
returns are and how long do we need to wait for the following close 
return (this is closely related to the gaps in the binary decomposition). 
The combination of these properties depends on quantitative properties 
of the fixed point. See Theorem I4.15[ 

The fact that the gaps in the binary decomposition play a role is 
clearly illustrated in the numerical calculations in [DEdlLOG] . 

4.6.2. The binary decomposition. Given a number n G N, we will use 
the binary expansion 

n = 2"'° + ■■■ + 2'"'^" (4.44) 

where mo > mi > . . . > m^^ > 0. 

If necessary, we will use mj{n) to emphasize the dependence on n of 
the power rrij in (14.441) . 

We will denote the integer part of any number u by [u]. Notice that 
that mo(n) = [log2(n)], and r„ < mo. 

Using (12.51) several times we have for all p > and for any x G [—1,1] 
that 
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A,(x,n) = J2 l^.,n(x)r A, [r-'^'--^'"''{x),2-^) (4.45) 



j=0 

When the initial value xq = x of x„ = satisfies |x| < |r.mo+i|, 
the points 

t;, = f'"^(t;,_i), v^i=x (4.47) 

where mo > . . . m^^ > 0, are return points to a small neighborhood of 
in the domain of renormalization. 

The decomposition in dyadic scales (14.451) will be our main tool for 
the analysis of the growth of (I4.45p . 

We will show in Lemma 14.131 that the size of the weights "^j^nix) 
are controlled using the renormalization. Then, using Corollary 14.81 we 
will get control over the size of the Ap(t;j_i, 2™^ ). Combining these 
results, the Lyapunov condition 12.91 will follow for the orbit of zero, 
thus establishing the central limit theorem. 

4.6.3. Estimating the weights. The main tasks of this section will be to 
estimate the weights in (14.451) and to show that they do not affect 
too much the conclusions of renormalization. 

For / G yVs{gk) fixed, recall the notation A/ = /(I), and for each 
nonnegative integer m, Tm = /2'"(0). Then 

> _ Tn+l 

At»/ - — — 

J- n 

If necessary, we will use r{ to emphasize the role of /. For the Feigen- 
baum fixed point g^ of order 2k, we use A^ = A^^. We have that 

n — -^k 

Recall that the map gk has a series expansion of the form 

oo 

gk{x) = l + Y,bix''' 

1=1 

and that g'k{l) = A^-^^ 

Lemma 4.12. Let hk be the function over the interval [—1, 1] defined 
by 

hkix) = ^ (4.48) 
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for X ^ and /ifc(O) = 2k bi. Then h is a negative even function, 
increasing on [0, 1]. 

Proof. From the fixed point equation 



9k{x) = -^gkigkihx)) 



we have that 



hkix) = '-^^^h{X,x) (4.49) 
9kW 

Repeated applications on (I4.49P gives 

h,{x) = h,iXr^x)fl'-^^^^ (4.50) 

Notice that hkiX^'^^x) — ^ hk(0) uniformly on [—1, 1] as ri — ^ oo. Fur- 
thermore, the factors in the product fl4.50p converge exponentially to 
1. Therefore, we can pass to the limit as n ^ oo in fl4.50p and get 



n=l 



Each term in the product is a nonnegative decreasing function on [0, 1] 
therefore, hj^ is increasing on [0, 1]. □ 

In the rest of this section, we assume that the domain of universality 
Ws{gk) is fixed. From the exponential convergence of T™/ to gk, we 
can choose a small neighborhood U of gk so that of 

T{Ws{gk) nU)c Wsigk) n U (4.51) 

The following result gives control of derivatives for points close to the 
origin. Some of the bounds for the derivatives in Lemma 14.131 are done 
in |VSK84j using the thermodynamic formalism. 

Lemma 4.13. Assume that f E Wsigk) and as before, let hk{x) = 
g'f,{x)x^~'^^ . For each finite decreasing sequence of nonnegative integers 
{rrijY-^Q, and \x\ < iFmo+il; let {vjYj^_^ with v^i = x be as in \4-47\ 
There exists < e ^ 1, depending only on U, such that if 

G = gk{Xk)-e (4.52) 

c = |MAfc)|-e (4.53) 

d = |/ifc(0)|+e (4.54) 

then 1 < cG"^^^^ < d, and 
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(1) For each < j < r, 

(2) Forl<j <r 



cG 



2fc-l 



mj_i 



2k-l 



< 



< d 



2k-l 



(4.55) 
(4.56) 



Proof. By (14.511) . we can assume without loss of generality that 

sup I A^m f — Xk\ < e 

\T'^f{Xk)-gk{Xk)\<e 
for some < e ^ 1. Since gk{Xk) > \Xk\ and 



(4.57) 
(4.58) 



> 1, 



by choosing < e ^ 1 smaller if necessary, we get 1 < cG"^^'^ < d. 

The exponential convergence implies that there is < 5 ^ 1 such 
that 



sup |r.m+ir„^| < |At'"/| < |Afc| + 5 



(4.59) 



We know that or all m G N, T"*/ is even, decreasing in [0,1], and 
sup^gj |T™/(t>)| < 1. Shrinking U if necessary, from (14.581) . we can 
assume without loss of generality that 

{9k{Xu)-t}<\T^f{x)\<l (4.60) 

for all x e [Afc — (5, |Afc| + 5]. 
(1) For j = we have that 

It follows from (1457]) that 

Then, by (14.601) we get 

{9k{Xk) - e}|rmo| < |t;o| < \^mo\ 

We continue by induction on j. Assume (14.551) holds for all numbers 
< ^ < J, where j < r. For j + 1 we have that 



^3+1 = f 



From fOTj) and fH39|) . it follows that 



{9k{Xk) -e}|r„^^J < \vj+i\ < |r^^._^j 
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(2) Shrinking U if necessary, we can assume that 



sup 



< e 



(4.61) 



for all f eWs{gk)nU. 

We know from Lemma [4.121 that 

MXk)\\x\''~' < \g'k{x)\ < \hm 

for all |x| < |Aa,.|. Hence, for each m G N, (I4.6ip and the exponential 
convergence of T"^f imply that 

(IMAfc)l -e) \x\''-' < |(r-/)'(x)| < (|/i,(0)| +e) [xp^-^ (4.62) 

forxG [/(Ay)|r„|,|r„|]. 

Therefore, from part (1) of the present Lemma, fl4.62p . and the iden- 
tity 

we have that 

2k-l 



/ 



< d\T-} v.^-^ < d 



r 



2k-l 



for all < j < r — 1. 



□ 



4.6.4. Lyapunov condition for all times for functions in the domain of 
universality. In this section, we show now that the Lyapunov condi- 
tion 12.91 holds along the whole sequence of positive integers for orbits 
starting at any point of the form x = /'(O), / G N, with / G Ws((y'fc) flW. 

Proposition 4.14. Let f G Ws(5'fc)nW be fixed, and for eachn G N, let 
mQ{n) = [log2(n)], r„ + l be the number of ones in the binary expansion 
ofn, and Imoin) be the interval [-\T,rno{n)+i\,\^mo(n)+i\]- 



(1) Let c, G, and d be as in Lemma 4-13. For any real number 
p > 0, there are constants ap, bp such that 
p 

Ap{x,n) 



< 



< br, 



P _)^(^2fc-l)mr„ / — ^^2^kp ^^-^rnoin) " ^ 1 ^{^2fc-l)m,., 
for allx G Imoin)- 

(2) For any real number p > 2, we have that 

Ap{x,n) 



(4.63) 



lim sup , 

"-°°xe/„o(„) {A2{x, n)Y'^ 







(4.64) 
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Proof. For each m G N, let = [— ir^+il, iTm+il]. Let n be fixed 
for the moment, and consider the sequence of powers, {mj}^'!^, in the 
binary expansion of n, see (14 .44^ . 

Let X G I-mo be fixed, and consider the sequence of returns {f 
with = X as in (I4.47p . 

From (I4.45P , we have for each p > that 

Ap{x,n) = J2 |^,,„(x)rAp(t;,_i,2'"0, (4-65) 

j=0 



(1) From Lemma [4.131 we have that 



2k-l 



where 1 < cG"^^ ^ < d. Since 

with < tUfe < 1, there is a constant a such that 

for all X G Imo{n)- By Corollary 14.81 there are positive constants Cp, Cp 
such that 



Cp < Ap(z/,„i,2" 



.p ^ ..pv-'.-i,- < Cp (4.67) 

Using (14.661) and (14.671) . we obtain estimates for the size of each term 
in (I4.65p . Since Xf^^Pp > 1, the estimate (I4.65P follows by noticing that 
mo — rrij > j. 



(2) Since < r„ < mo = [log2(r;,)], (I4.63P implies that there are con- 
stants < Cp < Dp such that for all n large enough 

CpiXfppr' < Apix, n) < Dp{rf%}"^« (4.68) 
From Proposition 14. 5[ we know that X'l'^^Pp > 1 for all p > 0, therefore 
lim inf Ap(x, n) = oo (4.69) 



Assume that p > 2. By Corollary 14. 8[ we have that for any given 
e > 0, there is an integer such that if m > M^, then 

Ap(a;,2")<e(A2(x,2"^))^'/2 



for all Ixl < |r^ 
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Denote by Sr„ = {0, . . . , r„} and let 

A^ i ^ Ap(M,2'"0 

= V " (A.(«.2".)K^ - ' 

and An = Sr„ \ A^. For each x G Imo and s = 2,p, we split As{x, n) as 

A,(x, n) = Hs{x, Al) + Al), 
where ^^), i = 1, 2 is defined by 

H,{x,J^n) = J2 |^,,n(x)rA.(t;,-i,2-0 

Notice that is the sum of bounded number of terms (at 

most Mg) all of which are bounded. Hence, 

Hp{x,Al)<a (4.70) 

for some constant > 0. 
In addition, notice that 

H,{x,Al) < e{H,{x,Al)r/' < e{A2{x,n)r/' (4.71) 

Combining iKTOh . fHTTj) . we get 

Ap{x,n) ^ Hp{x,Ai) + Hp{x,Al) 



(A2(x,n))p/2 - (A2(x,n))p/2 

< e+ ^' , 

Passing to the limit as n — > oo, we obtain from fl4.69p . that 

1- Ap(x,n) 
hmsup sup . . . ^ ^- 

By letting e ^ 0, we get (gSH). □ 

For any Z G N fixed and for all m > [log2(0]5 denote by Im = 
[—\Tm+i\, \Tm+i\]- We will show that the Lyapunov condition 12.91 holds 
for orbits starting at any point on {/'(0)}igN- 

Theorem 4.15. Let f G 'Ws{gk)^l^ be fixed, and I G N. For all p > 2, 

AJf(x),n) 

hm sup ^ \, = 4.72 

(A2(f(a;),n))^/^ 

The limit (I4.72p is not uniform in I. 
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Proof. The idea of the proof is to use the dyadic decomposition (I4.45P 
and show that for any fixed integer / > 1, the Lyapunov functionals 
Ap(/'(x), n) and Ap{x,n + I) for x G Imo{n+i) are very similar for large 
values of n (see f l4.77p ). 

Let n G N be a large integer and decompose n + / in its binary 
decomposition: 

n + l = + • ■ ■ 2™-- 

where mo = mo{n + /) = [log2(n + /)]. Consider x G lmo(n+i)+i, and 
observe that 

n+l 

j=i+i 

Using {rrij}, consider the sequence of returns {v}j with v_i = x (see 
(Hai])). It follows that 

r 

Ap{f{x),n) =\^n+i^o{x)\PAp{x, 2'"«) + J]|^„+,,(x)rAp(t;,_i,2"^0(4.73) 
where 

I 

A,(x,2™o) - A,(x,2™o) = J^lif"^""') ° /'(^)r (4-74) 



We will show that the difference fl4.74p is small compared to Ap[x, 2*"°). 
Indeed, since 

it follows from Lemma 14.131 that 



c 



U |(2fc-l)mo 



< 



(f-'-^)'°fW|< ,t.,„„ (4,75) 



where Cj and dj depend only on j. Since fl4.74p contains I terms, (14.171) 
and Corollary 14.81 imply that 

A,(x,2-°)-A,(x,2-°) 1 ^ 

Ap(x,2™o) - ''P^^2fcp^^)™o ^ ■ ' 

as mo oo. Therefore 

lim sup ^ = 1 (4.77) 

The limit (14721) follows from Proposition |411 and fli???]) . □ 
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We will use the result of Theorem 14.151 together with the results of 
Theorems 12 . 1 1 and 12 . 21 to prove the central limit theorem for systems in 
the domain of universality Ws{gk)- 



4.6.5. Proof of Theorem \2.^A By Theorem 12.11 and Theorem 14.151 it 
suffices to choose cr„ that satisfies fl2.10p . 

Recall that for any integer ra, mo{n) = [log2{n)] and r„ + 1 equals 
the number of ones in the binary expansion of n. 

For / G N fixed, denote by yi = /'(O). Let c, G, and d be as in 
Lemma 14.131 From (14.771) and Proposition 14.141 for each p > there 
there are constants Cp^i and Dp i such that for all for all x G Imo{n+i), 



(2fc-l)m, 



< 



n + l 



(Af%)™o(n+0 



< D 



Therefore, we have that 
{A{f'{x),n)r 

for some constant C > 0. 
Hence, if 



-1 



then, it suffices to consider 



with 7 > log2(p?/v^). 
In the case 



,Q2k 



d^Xf 



< 1 



> 1, 



it is enough to consider 



P,i 



I n + l 



x 



{2k-l)m, 



n + l 



Pi 



mo{n+l) 



(4.78) 



(4.79) 



(4.80) 



with 7* > hg.id^Xf'pl) + \og,{cG''^-'/^). 

To obtain the Berry-Esseen estimates of Theorem 12.21 we choose (T„ 
that satisfies (12.151) 

If f l4.79p holds, then it suffices to consider 



(Tr, 



1 



with e > log2(p?/v^) + log2(pi/P3) + 1- 
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In the case f l4.80p . it suffices to consider 

1 

With 9* > \og,{d^Xfpl) + log.icG^'^"'/^) + logM/p-s) +1 □ 

Remark 4.4. If we consider an increasing sequence of iterations {nk}k 
which have lacunar binary expansions, that is 

^"fc ^ g 

as — s> OO then, from (14. 78 p . it suffices to consider cr„j. = n'j^^^^^^ to 
have a central limit holds along {uk}. 

Remark 4.5. In the quadratic {k = 1) case, the power which arises 
from the argument in the central limit theorem presented here com- 
puted in |DEdlL06] . The result is 7 = 3.8836 . . . The method used in 
[DEdlLOG] follows a numerical scheme similar to [Lan82] to compute gi 
and the spectral radii of the operators JCg^^i and JCgj^^2- 

The numerics in |DEdlL06] give a con&mation that the dependence 
of the optimal cr„ on the number of iterations is a power of exponent 
7^, which is not too far from the the 7 coming from the arguments in 
this paper. 

4.6.6. Conjectures for orbits of points in the basin of attraction. The 
topological theory of interval maps |CE80l ICEL80j imphes that for 
each / G V n yVsidk), the set set Cf = {/"'(0)}„gf^ attracts all points in 
[—1,1], except for periodic points of period 2" - which are unstable - 
and their preimages. More formally, for any point a; G /, either 

(a) lim„^oorf(/"(a;),C/) = 0, or 

(b) /^"(/'(x)) = /'(x) for some nonnegative integers / and n. 

In the second alternative, we have | | > 1- 

The set C/ is commonly referred to as the Feigenbaum attr actor since 
it attracts a set of full measure (not open, however). It is also known 
that is uniquely ergodic since it is a substitution system [CEL80j . 
The set of points Bj that are attracted to Cf is called the basin of 
attraction of /. 

Numerical simulations |DEdlL06] for the quadratic Feigenbaum fixed 
point and the quartic map suggest that the effective noise of orbits 
starting in the basin of attraction Bf and affected by weak noise ap- 
proaches a Gaussian. We conjecture that the is indeed the case. 

Of course, the convergence to Gaussian is not expected to be uniform 
in the point in the basin. We have already shown that the convergence 
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to Gaussian is false for the preimages of unstable periodic orbits which 
are accumulation points for the basin. 

5. Central limit theorem for critical circle maps 

In this section, we consider another example of dynamical systems 
with a non trivial renormalization theory, namely circle maps with 
a critical point and golden mean rotation number. The t heory has 



been developed both heuristically and rigorously in |FKS82t IORSS83 
IMes84l [LiSM [SK871 [dF99] 

We will adapt the argument developed in section 14.11 to the case of 
critical circle maps with golden mean rotation number. 

5.1. Critical circle maps with golden mean rotation number. 

We will consider the following class of maps of the circle. 

Definition 5.1. The space of critical circle maps is defined as the set 
of analytic functions /, that are strictly increasing in in M and satisfy 

Gl. f{x + l) = f{x) + l 

G2. / has rotation number /? = ■ Recall that for a circle map 
/ the rotation number r(/) is defined by 

r{f) = hm 

n—ioo fl 

G3. (0) = for all < j < 2k, and /(2fc+i)(0) ^ 0. 

From the well known relation between the golden mean and the 
Fibonacci numbers, Qq = 0, Qi = 1 and Qn+i = Qn-i + Qn, it follows 
that QnP — Qn-i = (— 1)"~^/5" is the rotation number of the map 

/(n)(x) = /«"(a;)-Q„_i (5.1) 

The sequence of maps /(„) will be useful in section 15.31 where we 
analyze the cumulants of the noise. 

5.2. Renormalization theory of circle maps. There are different 
rigorous renormalization formalisms for circle maps. For our purposes, 
we will need very little about the renormalization group, so that we 
will use the very basic formalism of scaling limits. 

Definition 5.2. Given a map / as in Definition [5l let /(„) be as in (15.11) 

and denote A(„) = /(n)(0). The n-th renormalization of / is defined by 

fn{x) = fin) {\n~l)X) (5.2) 



58 O. DIAZ-ESPINOSA AND R. DE LA LLAVE 

Remark 5.1. Since the rotation number of f{n){x) is (— we 
have that 

for all n G N and a; G M. In particular, for x = we obtain that 
(— 1)"~^A(„) > 0. Therefore, it follows that each function fn{x), defined 
by fl5.2p . is increasing in x and satisfies fn{x) < x 

It is known |dFdM99l lYamOll IYam02l IYam03j that for every k E N 
there is a universal constant — 1 < < and universal function rjk 
such that: 

RCl. T]k is increasing, r]k{x) < x, and T]k{x) = Hi^{x'^^^^) for some 

analytic function Hk- 
RC2. The sequence of ratios a„ = A(„)/A(„-i) converges to a limit 

-1 < Afc < 

RC3. For some < 5^ < 1, \\fn{,x) - 'qk{x)\\co = 0(5^) where the 
norm || ■ \\co is taken over an appropriate complex open domain 
Dk that contains the real line M. 

RC4. The function rjk is a solution of the functional equations 

Vk{x) = jk^k ii^Vki^x)) (5.3) 

Vk{x) = ^Vk ihVkihx)) (5.4) 
RC5. The domain Dk can be taken so that dDk is smooth and 



XkDk C Dk (5.5) 
lim sup{| Im(z)| : z G Dk, \z\ > n} = (5.6) 



Xk'vkiXlDk) C Dk (5.7) 
XkVkiXkDk) C Dk (5.8) 

A few useful relations can be obtain directly from (15.31) and (15. 4p . 
For instance, by letting x = we get 

Vkil) = xl VkiXl) = xl 

Taking derivatives on ( 15. 3p and (15. 4p we have that 

V'kix) , / 1 .,2 ^ 

- Vk -^VkiXkX, 



v'kiXlx) " \Xk 

- VkiXkM 
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Then letting a: — > we get 

%(1) = 

Solutions of the equations fl5.3p . (15.41) are constructed in |Eps89] for 
all orders of tangency at the critical point. 

In the case of cubic circle maps, |Mes84l ILdlL84j using computer 
assisted proofs constructed a fixed point. The exponential convergence 
follows from the compactness of the derivative of the renormalization 
transformation at r]i, which is a consequence of the analyticity improv- 
ing of the auxiliary functions. This is part of the conclusions of the 
computer assisted proofs. 

5.3. Renormalization theory of the noise. In this section, we will 
develop in parallel two renormalization theories for the noise. For the 
purposes of these paper, either one will be enough so we will not show 
to which extent they are equivalent. We will assume that the order of 
tangency k is fixed. 

The renormalization group scheme for critical maps with rotation 
number (3 gives information at small scales and at Fibonacci times. 
Observe that if / is a circle maps, then Lyapunov functions A-^ and 
defined by (12. Sp and (12.20 are periodic with period 1. 

Denote by 

Since / is increasing, equation (12.50 for the Lyapunov functions Ap 
implies that 

knix) = (/(„_!) o /(„_2)(a;))^^n-2(a;) + A;„_i(/(n_2) (a;)) (5.9) 
knix) = (/(„_2) o /(„_i)(x))*'A;„_i(x) + A;„-2(/(n-i)(a;)) (5.10) 

5.3.1. Lindeberg-Lyapunov operators. To study the propagation of noise 
at small scales, we will use the following definition 

Definition 5.3. Let / be a critical circle map, and for each p > and 
n G N, define the operators 

Un,ph{z) = [f^__2{an-2fn-l{0in-lz))Y h{an-iz) 
Tn,pq{z) = [f^_^[a~^2fn~2{(yn~ian~2z))Yq{an-l(yn~2z) 
Rnh{z) = /l(a^l2/n-2(«n-l«n-22)) 

PnQiz) = g(a„_2/„-i(a„-i2;)) 
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acting on the space or real analytic functions in the domain D^. For 
p > 0, and n E N, the Lindeberg-Lyapunov operators, ]Cn,p and )Cn,p, 
acting on pairs of real analytic functions in Dk are defined by the 
matrices of operators 

l^n,p — I J g I r^n,p — I j g I l^-^'-'-J 

where / is the identity map and is the zero operator. Similarly, 
consider the operators acting on the space of real analytic functions on 
defined by 

Uph{z) = [n'f.iXkVkhz))]^ h{\kz) 

Tpq{z) = [vUK'Vk{Xlz))Y qiXlz) 

Rh{z) = hiX.^iXlz)) 

Pq{z) = q {XkT]k{Xkz)) 

The Lindeberg-Lyapunov operators ICp and JCp are defined by 

^p={f I') = ( ^ ) (5-12) 



Notice from ( 15. 9p and ( IS.lOp that the change variables 

X(n-l)Z, kn{z) = kn{X^n-l)Z), 



X 



implies that 



k 

kn-1 
k 

kfi-l 



^n,p ' 
^n,p ' 



3,p 



3,p 



k2 

~ki 
~ki 



(5.13) 
(5.14) 



Remark 5.2. Notice that the derivatives involved in Definition 15.31 are 
positive. Therefore, the notions of cumulant operators and Lindeberg- 
Lyapunov operators coincide in this case. 

Equations fl5.13p and fl5.14p are the analogs to (14. 6 p for period dou- 
bling. These equations measure the growth of the linearized propa- 
gation of noise at Fibonacci times for an orbit starting at the critical 
point 0. 
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5.3.2. Exponential convergence of the Lindeberg-Lyapunov operators. 
An important consequence of the exponential convergence of to rj^ 
is that the Lindeberg-Lyapunov operators )Cn,p and /C„^p will converge 
exponentially fast to operators JCp and /Cp respectively, as n —>■ oo. 
Indeed, we have that 

Lemma 5.1. Let f be a circle map of order k satisfying G1-G3. 
The operators /C„,p, ICp, /C„,p and fCp acting on the space of pairs of 
bounded analytic functions defined on some compact set C that 
has smooth boundary and contains [—1,1] in its interior, are compact. 
Furthermore, for p fixed, 

- VkWsk + Wfn-i - Vkllsk) (5-15) 

\\K:n,p - iCpW < CpiWfn - r]k\\Bk + Wfn-l - VkWBt,) (5.16) 

Proof. We only show (15.161) . since the proof of fl5.15p is very similar. 

For any A G C, recall the notation A"^ = {z E C : d{z,A) < e}. 
Since — \k\ = 0{S'^) for some < 5^ < 1, we have that for any 

compact Bk C Dk and for all n large enough 

Oin-iBk C XkBk C Dk 

where e > is small. 

For z G Dk, we have that 

\(^n—2fn— 

+ \an-2 - \k\\r]k{an-iz)\ + \Xk\\r]k{an^iz) - r]k{\kz)\ 

Therefore, using Cauchy estimates, we have that for any compact set 
Bk d Dk containing 

||a„_2/„-i o an-i - XkVk o Afellsfc < C (||/„-i - Vkllsk 

+ \\fn-2~Vk\\Bj (5.17) 

Therefore, we by taking e small enough, we have that for all n large 
enough 

a„_2/„-i(a„-i5fc)' C XkVkihBkf C Dk (5.18) 

Consider Bk with smooth boundary and let Wk = Xki^ki^kBk) ■ For 
z E Bk we have 

l/n-2("n-2/n-l(an-l^)) - r]'^iXkVkXkZ))\ < | /^_2 (a„-2/„-l («n-l2;) ) 

fn-l{an~iz))\ 

+ \r]'k{an-2fn-i{an-iz)) - if]'k{\kr]kXkz))\ 
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Then, from flS.lSp and fl5.17p we get 

ll/n-2 O («n-2/n-l O On-l) - % O (Afe// O Afc)|| < |/^_2 - V'kWu 

+ C'hfclkJI (ll/n-1 - VkWsk + ll/n-2 - ^fclsj 

Using Cauchy estimates we obtain 

ll/n-2 ° (an-2/n-l O On-l) - % O (AfcT/ O Afc)||iJ, < C - T^fcHz), 

+ ||/n-2-%|Dj (5.19) 

Denote hy O = 'C \ {x + iy : x < 0, ?/ = 0}, and let C -D^ be a 
bounded domain containing [—1, 1] with smooth boundary such that 



From f l5.19p we can assume, by taking e smaller if necessary, that 



an-2fn^l{,an-lBk) C Xk^kiXkBh) C O 

Recall the operators t/„ p. Up, Pn,p and Pp defined in Definition 15. 3[ We 
restrict these operators to the space of bounded analytic functions on 
Bk. Denote by = XkVkiXkBk) and let ip^z) = defined on F^. 
Hence, using the triangle inequality and Cauchy estimates, we get for 
any bounded analytic functions h, q on B^ that 

\\{Un,p - Up)h\\B, < ||^||yJ|/l||A,i?J«n-l - Afc| 

+ C'||V^'||yJ|/i||A,B, (II/.-1 - Vkh, + ll/n-2 - Vk\U) (5.20) 

\\{Pn,p-PMBk<C\\q'\\wA\\fn~l-r]k\\B, 

+ \\fn-2-Vk\\Bk (5-21) 

Combining (15:201) and (ICTD we obtain (15J6|) . 
For zi, Z2 & Bk close enough we have that 



Tn,ph{zi) - Tn,ph{z2)\ < [ip O f'^_^{an-2fn-l ° "n-l))' 



'-1 — ^2 



k 



\Pn,ph{zi) - Pn,ph{z2)\ < 1 1 (g o (a„_2/„-i o a„_i ) )' 1 1 | Zi - 2:2 1 

Therefore, the compactness of the operators ICn,p follows from Cauchy 
estimates and the Arzela-Ascoli theorem. □ 



As in section 14.21 the spectral properties of the operators /C„ p, ICp, 
fCn,p and Kp will prove important for our analysis of the Gaussian 
properties of the scaling limit. 

Notice that for p G N, since fn and rjk are increasing, the cumulant 
functions defined by fl3.43p and the Lyapunov functionals (12.20 coincide. 
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Hence, the study of the properties of cumulants is equivalent to that 
of the properties of the Lindeberg-Lyapunov functionals. 

5.4. Spectral analysis of Lindeberg-Lyapunov operators. From 
Lemma [5.11 we know that the Lindeberg-Lyapunov operators /Cp, /Cp 
when defined on in a space space of real valued analytic functions with 
domains in D that satisfy ( 15. 5p . (15. 6p . (15. 7p . ( 15.80 are such that their 
squares are compact. Again, the compactness follows from the fact 
that the operators f/„.p, r„ p, Q„ in Definition 15.31 are analyticity 
improving. 

The compactness of the operators /C„,p, /C„,p, is slightly more subtle 
since they have a block (in the lower diagonal) that is the identity. Us- 
ing the fact that this block appear in the lower diagonal, the operator 
can be made compact by considering that the second component is an- 
alytic in a slightly smaller domain. In this way, the identity becomes an 
immersion, which is compact. If the domain for the second component 
is chosen only slightly smaller, the analyticity improving properties of 
the other operators are still maintained. 

Remark 5.3. An slightly different approach is to remark that the 
square of the operators consist only of compact operators. For our 
purposes, the study of the squares will be very similar. 

Notice that the operators /Cn,p, ^n,p preserve the cone C of pairs of 
functions, both of whose components are strictly positive when re- 
stricted to the reals. Note that the interior of the cone C are the pairs 
of functions whose components are strictly positive when restricted to 
the reals. 

It is also easy to verify the properties (i), (ii), (iii) in Section 
Again, (i), (ii) are obvious and we note that if a pair is not identically 
zero, then, one of the components has to be strictly positive in an 
interval. The structure of the square of the operator implies that the 
component is strictly positive in an interval which is larger. This, 
in turn, implies that the second component becomes positive. Both 
components become strictly positive by repeating the operation a finite 
number of times. 

Therefore, we can apply the the Krem-Rutman theorem and obtain 
a result similar to Proposition 14.21 

Proposition 5.2. Let ICn,p and ICn,p be the Lindeberg-Lyapunov oper- 
ators defined by ( 15. lip . Denote by /Coo.p = ^'Cp and /Coo,p = ^p- For all 
p > and n E NU {oo}, let pn^p (pn,p) be the spectral radius of ICn^p 
(K,n,p)- Then, 
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0^) Pn,p (Pn,p) is a positive eigenvalue of lCn,p (^n,p)- 

b) The rest o/Spec(/C„,p) \{0} (Spec{ICn,p)\{0} ) consists of eigen- 
values fi with < Pn,p (\p\ < Pn,p)- 

c) A pair of positive functions {ipn,p, (pn,p) is an eigenvector of }Cn,p 
(JCn.p) if and only if the corresponding eigenvalue is Pn.p (f^n,p)- 

5.4.1. Properties of the spectral radius of the Lindeberg-Lyapunov op- 
erators. The convexity properties of the Lyapunov functions Ap will 
also imply similar properties for the cumulants ICn,p and /Cn,p- In the 
next Theorem, we will only consider the the operators )Cp and Kp. We 
also remark that the argument in Proposition 14.41 applies in our case 
also, so that the spectrum is independent of the domain considered. 
Hence, we will not include the domain in the notation. 

Theorem 5.3. Let pp be the spectral radius of ICp. ( Similarly for ^p 
and JCp.) 

Kl. For allp>0 

^f'Pp > 1 

K2. The map p ^ pp is strictly log-convex 

K3. The map p i— > \og{pp)/p is strictly decreasing. 

A Similar result holds for pp. 

Proof. (Kl) For any pair of positive functions [h, £], we have that 

ICf'ih, I] > [T^h,T;'i] 
If [h, i] is the dominant eigenvector of /Cp, then we have that 
pf^[h{z), i{z)] > [T™/i,T™£](z) 

m 

Since ?7^(1) = A^^^, the conclusion of (Ml) follows by letting z = 0. 
To prove [K2.] [K3.] we proceed in a similar way as in the proof of The- 
orem |431 We study the eigenvalue equation for the leading eigenvalue 
and use it to construct positive functions that satisfy eigenvalue equa- 
tions with a positive remainder and then, we apply Proposition 14.61 
We suppress the subindex n for typographical clarity. We will use the 
notations for operators introduced in Definition 15.31 

The eigenvalue equation for the leading eigenvalue of /Cp is 

Ppip = Rhp + Tp£p (5.22) 
Ppip = hp (5.23) 
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We raise the equations (15.221) to the power a > 1 and use the inequal- 
ity for the binomial theorem and the elementary identities {Tpip)"" = 
Tpa^p, {Rhp)°' = Rhp. Then, we have: 

p^rp = {Rhp + Tpipr (5.24) 

> {Rhpr + {Tpipr (5.25) 

= Rh'^ + TpJ^ (5.26) 

Pp^p = h; (5.27) 



Hence, we obtain 



and we have that the inequality is strict in the first component. There- 
fore 

}cip[h^,n>pl-[h-,n 

Using that the leading eigenvalue of JCp is the square of the leading 
eigenvalue of /Cp, we obtain, applying Proposition 14.61 that 

Ppa < Pp 

which is property [K3.]. 

To prove property [K2] , we multiply the eigenvalue equation for lead- 
ing eigenvectors of /Cp, ICp and raise to the 1/2 power. We use the 
identities {Rhp){Rhg) = R{hphg), [R{hphq)Y/^ = R [{hphg)^/^] as well 
as similar identities for the operator T. 

{pppgy^^hph.y/^ = [{Rhp + Tpip) ■ {Rh, + T/,)] 
> {RhpRhp + TpipTgigY^^ 

Multiplying the second components we have 

{KKf/' = {ppp.y^'iw^' 

In other words, we have 

/C(p+,)/2[(^^)^/^(v.)'/'] < m'^'ukky^mw^'] 

with the inequality being strict in the first component. Proceeding as 
before, we obtain 

P{p+q)/2 < {ppPqf^ 

which is equivalent to strict log-convexity. 

□ 
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5.5. Asymptotic properties of the renormalization. A conse- 
quence of Lemma EU] is that the asymptotic properties of the operators 
fCn,p and fCn,p become similar to those of Kp and Kp respectively. 

The following result, Corollary 15.41 is obtained as an application of 
Proposition 14. 7[ Since the proof is very similar to that of Corollary 14.81 
for period doubling, we will omit the details. 

Corollary 5.4. Let lCn,p, ^p, O'nd Pp be as in Proposition \4 ■ S\ Then 

1. There is a constant Cp > such that for all positive analytic 
pairs of functions [h, q] in the domain V of the Lindeberg- 
Lyapunov operators 

c;'p;<ICn,p---IC,,p[hiz),q{z)]<Cpp; (5.28) 

2. For p > 2, we have pp < (p2)^^^ 

3. For any p > 2 

A/(x,Q„) 
lim sup 7 = U 

™.e/„{A^(a;,Q„)}p/2 

where Im = [-|A(„)|, |A(„)|]. 

5.6. Proof of the central limit theorem for circle maps (Theo- 
rems 12.41 and 12.51) . In this section, we will prove a central limit for 
critical circle maps (Theorem 12.41) . The method of the proof is similar 
to the method developed in Section 14.61 for period doubling. That is, 
we will verify the Lyapunov condition (12.91) along the whole sequence 
of integers and then use Theorems 12.11 and 12. 2[ 

The renormalization theory developed in Section [53] gives control on 
As(a;, Qm) for all x sufficiently small (|x| < |A(m+2)|), see Corollary 15.41 
The main tool to obtain control of the noise on a segment of the orbit 
will be to decompose the segment into pieces where renormalization 
applies, that is segments of Fibonacci length. 

This decomposition is develop in Section 15.6. 1[ The main conclusion 
of this section is that the effect of the noise over a period of time n 
equals the sum of Fibonacci blocks - where renormalization applies - 
with some weights. See (15.301) . This weights measure how the effect of 
Fibonacci blocks propagates to the end of the interval. 

The effect of weights on the Fibonacci blocks is measure in Section 

5.6.1. Fibonacci decomposition. Given n G N, it admits a unique de- 
composition 

n = Qmo + --- + Qm.„ (5.29) 
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where mo > . . . > m.,.^ > are non-consecutive integers and Qm is the 
m_,— th Fibonacci number. If necessary, we will use mj{n) to emphasize 
the dependence on n. Notice that r„ < mo{n) < [log^-i n + 1]. 
We have the following Fibonacci decomposition for Ap{x, n) 

A,(x,n) = |vl>„„(x)r Ap (^"^'"^"■■■"'^'"'■(a;),Q„^,) (5.30) 

where 

= (/"-Q'"o--Q-.)'o (/«-.+-«'"o) (a;) (5.31) 

Since is periodic of period 1 for any m G N, Ap(x, n) can be 

expressed in terms of the functions /(„) defined by ( 15.11) . Indeed, let 

= f{mj){Vj-l) V-i = X (5.32) 

then, 

l=j+l 



and 



Ap(x, n) = Y,{'^n,,{x)YK{'^j~i^ Qm,) (5.33) 
i=o 

Let us fix a critical map of the circle of order fc, and let < e <^ 1. 
For all n large enough we have that 

ll/n - ^fcllc(Dfe) < e (5.34) 
Using renormalization, we will control the size of the size of the weights 

'^j,n{x) for \x\ < |A(,„o+2)|. 

5.6.2. Estimation of the weights. In this section we estimate the weights 
'^n,j defined in fl5.33p . The method used is very similar to to one de- 
velop for period doubling. 

Recall that r]k{x) = + Ylf=i^j^'''^^~^''^ ^ With hi > 0. Consider the 
function defined by 



and define 

Sk = inf Hk{x) Uk = sup Hk{x) (5.35) 

{x:|x|<A|} {x:\x\<\l} 

It follows that 

cx^^ < v'kix) < dx^^ (5.36) 
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In the rest of this section, we will assume with no loss of generality 
that all n are large enough so that 



sup \fn{x) - rjk{x)\ < e 

{x:\x\<l} 



(5.37) 



with < e < 1. 



Lemma 5.5. Let {rrijYj^Q be a decreasing sequence of non-consecutive 
positive integers. For \x\ < X{mo+2), l^t {vjYj=_i be defined by (15.321) . 

(1) For all < j < r for which rrij is large enough 



|A|||A(^^_i)| < It;, I < \Vki-\l)\\\m,~i)\ 



(2) Let c = Sk — € and d = Uk + e. Then, 

2k 



cX 



6k 



A(mj_i-1) 



< (/H,))'(t^,-i) < d 



Proof. (1) For j = 0, we have that 



A(m,_i-1) 



A 



(m,-l) 



2k 



V-1 



(mo-l) 



The exponential convergence of fm to rjk implies that 



V-l 



< 



A 



(nio+2) 



A 



(mo-l) 



< |Afcp + e 



A (mo-l) 

for some < e ^ 1. Hence, 

(|Afe|'-e)|A(„„_i)| < |t;o| < \r]k{-Xl)\\Ximo-i)\ 
By induction, assume that 

\vj-i\ < \vk{-Xl)\ |A(„^_,_i)| 



Notice that 



f{mj){Vj-l) — A(mj_i)/„ 



Since m,_i — rrij > 2, we have that 



X 



(nij-l) 



<M-Xl)\iXl + e)<Xl 



since e -C 1. Therefore, we have that 



(5.3^ 



(5.39) 



(5.40) 



|Afc|^ |A(mj_l)| < \Vj\ < IVki — Xl)] |A(mj_l)| 
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(2) Taking the neighborhood of the fixed point smaller if necessary, 
assume that for all m large enough 



sup 



X 



2k 



X 



2k 



< e 



Notice that 



A 



(mj-l) 



If TTLj is large enough, by the exponential convergence of the renormal- 
ized functions to the fixed point, we have that 



Therefore, from (15.361) we get 



cX 



6k 



\ \ 2fc 

A(mj_i-1) 

A(mj-l) 



< |(/(n^,)y(^.-l)| <^ 



\(m 



-1) 



A(mj-l) 



2k 



□ 



5.6.3. Lyapunov condition for critical circle maps. In this section we 
show that the Lyapunov condition 12.91 holds for subsequences of the 
orbits starting at points of the form x = /'(O), / G N. The main re- 
sult of this section is Theorem 15.71 which proves that the Lyapunov 
condition holds always if we choose appropriately subsequences (they 
are subsequences of numbers with very few terms in the Fibonnaci ex- 
pansion. The precise conditions depend on numerical properties of the 
fixed point. We show that if the fixed point satisfies some conditions, 
we can obtain the limit along the full sequence. 

First, we use Lemma 15.51 to estimate the growth of the Lyapunov 
functions Ap(0,'n,) at zero at all times. 



Proposition 5.6. Let c and d as in Proposition \5.5l For each integer 
n, let Qmo{n) be the largest Fibonacci number in the Fibonacci expansion 
ofn and letrn be the number of terms in expansion (I5.29P ofn. Denote 

by Imo = [— |A(mo+2)|, |A(mo+2)|]- 

1. For each p > 0, there are constants ap and bp such that 



y2km.rn 



< 



< br, 



d'" 



X 



(5.41) 



for all X G /, 



mo(n) ■ 
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2. If {rii} is an increasing sequence of integers such that 



lim , = (5.42) 

i^oo mo[ni) 



then 



V A(m,_l) J 



lim sup J^f^^l^^ = (5.43) 

J^°°^e/™o(",) (A2(a;,nj-))^/^ 

3. // the following condition 

{skXfYxf'pp>l, (5.44) 
holds (sk defined as in (15.351) ) then, for each p > 2 

AJx,n) 

lim sup \/ = 5.45 

(A2(a:,n))^/2 ^ ^ 

Proof. For a given n G N, let {mjYj'^Q be the decreasing sequence of 
numbers that appear in the Fibonacci expansion (I5.29p of n. 
(1) Recall the weights \l/„,j defined in (15.331) . By (I5.39P in Lemma [5751 
and (15.371) . we can assume without loss of generality that all mj are 
large enough, so that 

2k 

I < d"-' 

\(mj-l) 
for all X e Imoin)- 

The exponential convergence of /„ to rjk, implies that there is a 
constant C > such that 

C'\cXfY-'^^^ < Af ^"''■""'"^^^„j(x) < Cd^--^+^ (5.46) 

By Corollary 15.41 there is a constant D > such that 

< p;"^ Ap(t;,_i, Qm,) < D (5.47) 

for all j = 0, . . . r„. 

Since each term in the decomposition (I5.33P of Ap(0,n) is positive, 
for some constants and bp we have that 

(•V mo \ 2fcp Jn y ,nij \ 2kp 

for all X G Imo{n)- We obtain 05.411) from fl5.46p . 05.471) by noticing that 
niQ — rrij > 2j . 

Notice that if {ui} is a sequence of integers that satisfy (I5.42p then 
lim inf Ap{x, Ui) = oo (5.48) 
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On the other hand, if the condition (15 .44^ holds, then {c\^Y)^^^ pp > 
1. Therefore, 

hm inf Ap{x,n) = oo (5.49) 

n->oo xe/„Q(„) 



(2) Corollary 15.41 imphes that for any e > 0, there is an integer such 
that if m > M^, then 

Ap(a;,Q™) < e (A2(x, g„)f 
for all X G [— |A(m)|, |A(m)|]- Let Sr^ = {0, . . . ,r„} and define 

Ai = I J G : sup ^'^'^'"-^ < e 

I (A2(x,g^,))"/' 

and An = \ For each x G Imo and s = 2,p, we split As{x, n) as 

A,(x, = /7,(x, Al,) + if, (a;, ^^), 

where Hs{x, A\^, z = 1, 2 is defined by 

H,{x,A^^) = |^,,„(x)|"A,(t;,_i,Q^J 

Notice that iip(x, is the sum of at most bounded terms. 
Therefore, for some constant > 

Hp{x,Al)<C, (5.50) 

For the term Hp{x,A}ri) have 

Hp{x,Al) < e{H2{x,Al)r/' < 6(A2(x, n))^/^ (5.51) 

Therefore, combining (15.501) . (I5.5ip we get 

Ap{x,n) ^ C, ^ ^ 



(A2(x,n))p/2 - (A2(x,n))p/2 
Limit fl5:i3|) follows from fICTD . 

(3) If condition flSlliD holds, then limit (ISllSj) follows from ([539]). □ 

Remark 5.4. One natural hypothesis that implies condition (15.441) 
holds is that Sk = Kk{\l). In this case, H{\l)\f = 1. Then, (ICTD 
follows from Proposition 15.31 
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Since 



for all m, it follows that mo{n) ~ [log^-i(?2)] + 1 as — >■ oo. 

In the following result, we show that the Lyapunov condition [2l 
holds for all orbits starting at points of the form x = /'(O), / G N. 

Theorem 5.7. 1. // {rii} is an increasing sequence of integers 
that satisfy fl5.42p . then 

Ar,( f'(x), nj) 

lim sup ^^■^ ^ ^' ^ = 5.52 

^^o^.ei^,,.^^,iA,if^ix),n,)f' 

for all p > 2. 
2. Under condition (15 .44^ 

lim sup Mf\-),n.) (553^ 
for all p > 2. 

Proof. We will use a method similar to the one used in Theorem 14.151 
since the proof of (1) and (2) are very similar, we will only prove (2). 
Let Z e N be fixed. For any n G N, consider 

ri -\- I QniQ Qnir 

where {rrij} is as in (15.291) . For x G Imo{n+i), define the sequence of 
returns {vj}j with t;„i = x as in (I5.32p . Notice that 

Ap(/(x),n) = \^n+i,oi^)\''Ap{x,Q^,) + 

r 

Y.l'^n+lA^WM^i-uQmJ (5.54) 
i=l 

where 

I 

Ap(x,Q„o) - Ap(x,Q„J = J2 I (/'^'""~') ° /'Wr (5-55) 

i=i 

From the identity 



fP"o-.YoP(x) = i^^:^ 



and Lemma [5.51 it follows that 



Cj 



(/'3'"o-.yo/^(x)|<-A_ (5.56) 
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where Cj and dj depend only on j. Since ( 15351) contains / terms, Theo- 
rem (I5.3p and Corollary 15.41 imply that 

Ap(x, Qmo) ^p{x, Qmo) ^ Q ^ ^ g 

Therefore, 

lim sup — — = 1 (5.57) 

and the limit (15.531) follows from Proposition (15.61) . □ 

5.6.4. Proof of Theorems \2.4\ and \2.5[ Recall that for large integers n, 
mo{n) ~ [log^-i(?T,)] + 1 and r„ + 1 equals the number of terms in the 
Fibonacci expansion of n. 

By Theorem 12.11 and Theorem 15. 7[ it suffices to choose cr„ that sat- 
isfies f i230|) . 

Let c and d be as in Lemma 15.51 From (I5.57P and Proposition 15. 6[ 
for each p > there there are constants Cp^i and Dp i such that for all 

for all X e Imo(n+l), 

Therefore, we have that 

for some constant C > 0. 

If we consider an increasing sequence of iterations {nk}k which have 
lacunar Fibonacci expansions, see (15.421) . Then, from (I5.58P it follows 
that it suffices to consider (T„^ = ^'^^'^^ to have a central limit holds 
along {rifc}. 

Otherwise, under condition (I5.44p we have the following cases: 
Case (a) If 

d' 



< 1 (5.59) 



then, it suffices to consider 



with 5 > log^-i(p?/^). 
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Case (b) If 



it is enough to consider 



>1, 



(5.60) 



n' 



<5*+l 



with 6* > log^-i(rf3p?) + log^-i(cAf /yp^). 



The Berry-Esseen estimates follow by choosing (T„ that satisfies fl2.15p 
in Theorem 12. 2[ 

If f l5.59p holds, then it suffices to consider 



1 



with r > logp i {pi /y p^) + log^-i (p^/ps) + 1. 
In the case f l5.60p . it suffices to consider 



CTr. 



n' 



with T* > log^-i(d3p3) + log2(cAf /yp^) + log2(pi/p3 



-1 



□ 



Remark 5.5. We know from the theory of critical circle maps that 
{/'(O) mod 1 : / G N} is dense in T^. Numerics in |DEdlL06] suggests 
that the Lyapunov condition holds for all points in the circle. We con- 
jecture that this is indeed the case and that that the effective noise of 
orbits starting in the any arbitrary point of approaches a Gaussian. 
We expect the the speed of convergence in the central limit theorem to 
be not uniform. 



6. Possible extensions of the results 

In this section, we suggest several extensions of the results in this 
paper that are presumably be accessible. 

1. Assume that is a sequence of independent random variables 
with mean zero p finite moments such that 

with some a± in a small range. 

2. Assume that the random variables are weakly correlated 
(e.g. Martingale approximations). 
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These assumptions are natural in dynamical systems applications when 
the noise is generated by a discrete process. That is 

Xn+l = f{Xn) + (y^iVn) 

Vn+i = h{yn) 

is an expanding map or an Anosov system. 

Related to the central limit theorem (even in the case indepen- 
dent random variables (^„) of comparable sizes), it also would 
be desirable to obtain higher order asymptotic expansions in 
the convergence to Gaussian, namely Edgeworth expansions. 
We note that the estimates for the asymptotic growth of the 
variance of the effective noise ( (14.231) with p = 2) for systems at 
the accumulation of period doubling are obtained in |VSK84j 
using the Thermodynamic formalism. We think that it would 
be very interesting to develop analogues to the log-convexity 
properties of the Lindeberg-Lyapunov operators or the Edge- 
worth expansions with the thermodynamic formalism. 
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